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MODELING  MODAL  TRANSIENT  BEHAVIOR  OF 
GaAs  LASER  DIODES  WITH  TEMPERATURE 
DEPENDENT  RATE  EQUATIONS  AND 
COMPARISON  TO  EXPERIMENTAL  DATA 

CHAPTER  1 
INTRODUCTION 

Recent  advances  in  GaAlAs  laser  diodes  have  led  to  the  development  of 
low  cost,  long-lived,  high  power  devices  which  have  found  application  as 
transmitters  in  various  communication  systems.^  One  problem  with  the  use  of 
gain-guided  laser  diodes  in  communication  systems,  however,  arises  because 

significant  power  sharing  between  the  axial  lasing  modes  can  occur  over  a 

period  of  order  100  ns  when  the  lasers  are  pulsed.  These  modal  "instabilities" 
(known  as  mode  partition  noise)  can  cause  distortions  in  an  optical 
communication  system  if  the  receiver  samples  only  part  of  a  dispersed  optical 
beam.  The  power  in  the  different  axial  modes  will  shift  with  time  during  the 
pulse,  and  so  will  the  received  signal  if  only  selective  modes  are  sampled. 

This  problem  was  investigated  by  first  taking  experimental  data  on 
several  types  of  laser  diode  structures  to  quantify  the  magnitude  of  the  modal 
instabilities.  The  modes  are  dispersed  with  an  cchelle  grating  and  examined 
individually  with  an  avalanche  photodiode  (APD)  under  various  modulation 

conditions.  The  data  was  then  compared  to  predictions  of  the  modal 

instabilities  made  by  solving  temperature-dependent  rate  equations.  The 
thermal  impedance  of  the  laser  diodes  and  the  temperature  changes  in  the 
active  area  during  a  current  pulse  were  modeled  with  a  finite  element  code.  A 
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detailed  analysis  of  the  laser  was  done  2-dimensionally,  with  a  3-dimensional 
analysis  performed  to  determine  the  thermal  impedance  of  the  heat  sink.  The 
thermal  properties  were  then  entered  into  a  program  which  numerically 
solves  the  rate  equations  of  the  diodes  so  that  the  power  from  the  individual 
modes  of  the  diode  are  determined  as  a  function  of  time  into  the  current  pulse. 

Broad-stripe,  narrow-stripe  and  10-stripe  laser  diodes  were  examined  in 
this  thesis.  The  broad-stripe  and  10-stripe  laser  diodes  are  high  power  lasers, 
producing  200  to  600  mW  of  power.  The  10-stripc  laser  diode  array,  consisting 
of  10  narrow  stripe  lasers  spaced  close  together  on  the  same  chip,  is  an  attempt 
to  retain  the  high  power  of  the  broad  stripe  laser  while  reducing  filamenting, 
a  condition  where  lasing  does  not  occur  evenly  throughout  the  gain  area  but 
tends  to  concentrate  in  certain  areas  called  filaments.  The  narrow-stripe 
lasers  are  lower  power  but  were  also  examined  in  this  report,  mainly  as  a 
comparison  to  the  other  structures  and  to  help  better  understand  laser  diode 
behavior  in  general. 

This  report  is  organized  with  the  experimental  data  described  first. 

Then  details  are  given  on  the  finite  element  thermal  modeling.  The  rate 
equation  model  is  then  described,  and  the  lasers  are  modeled  with  the  rate 
equation  model  first  without  and  then  with  thermal  effects.  Finally  a 
comparison  is  made  of  the  experimental  and  calculated  results,  and  suggestions 
are  made  for  optimum  laser  designs  that  will  minimize  the  modal  instabilities 
of  the  laser  diodes  during  modulation. 
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CHAPTER  2 

EXPERIMENTAL  MEASUREMENTS  CHARACTERIZING  THE 
PULSE  RESPONSE  OF  LASER  DIODES 


2.1  Determination  of  the  pulse  response  of  the  modes  of  the  laser 
diodes  with  an  APD 

The  pulse  response  of  several  laser  diodes  were  examined  with  an  APD. 
Figure  2.1  shows  the  experimental  set-up  used.  The  laser  diode  output  is 
focused  onto  a  slit  at  the  spectrometer  entrance.  The  incident  beam  then 
illuminates  an  echelle  grating.  The  grating  separates  the  laser  modes  which 
can  be  either  imaged  with  a  CCD  camera  or  projected  onto  a  scanning  APD. 


The  APD  output  is  amplified  and  then  sent  to  a  Boxcar  integrator  which 
averages  the  data  over  a  desired  time  period. 


FIGURE  2.1  Experimental  setup  used  to  measure  laser  diode  pulse  response. 
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Several  types  of  lasers  were  examined  in  this  report.  Table  2.1  lists  the 


various  lasers  along  with  their  pertinent  characteristics. 

TABLE  2.1 

Lasers  Tested  Experimentally 


Designation 

Structure 

Type 

Width 

Threshold 

pm 

pm 

m  A 

SOW-5-300 

5 

300 

75 

SOW-60-300 

broad  stripe 

300 

97 

SOW-5-600 

narrow  stripe 

5 

600 

130 

SOW-60-600 

60 

600 

102 

MOW-4-250 

4 

250 

50 

MOW-100-250 

■ESSSSHI 

lO-stripe 

100 

250 

221 

The  first  four  lasers  in  Table  2.1  are  AlGaAs  graded-index 
separate-confinement  heterostructure  (GRINSCH)  single-quantum-well  (SQW) 
lasers  mounted  p-side-down  from  McDonnell  Douglas  Electronics  Systems 
Company. 2.3  These  lasers  are  gain-guided  in  the  lateral  dimension.  The 
10-stripe  multi-quantum-well  (MQW)  laser  is  from  Spectra  Diode  Labs  (SDL) 
and  has  4  p.m  stripes  on  10  pm  centers.'*  The  narrow  stripe  MQW  laser  is  also 
from  SDL  and  has  a  4  pm  wide  stripe.  Figures  2.2  and  2.3  show  the  10-stripe 
laser  diode  array  and  the  broad  stripe  laser  diode  structures.  The  wavelengths 
are  all  near  820  nm. 

Figure  2.4  shows  the  output  waveforms  for  twelve  dominant  modes  of 
the  SQW-60-300  laser  modulated  at  105+100  mA  with  an  80  ns  pulse  with  a  50% 
duty  cycle.  At  this  modulation  depth  the  laser  produces  80  mV/  of  light.  For  the 
broad  stripe  lasers  another  cylindrical  lens  was  used  to  focus  the  broad  stripe 
onto  the  detector. 

Two  types  of  phenomena  cause  modal  instabilities  in  pulsed  laser  diodes; 
mode  buildup  transients  and  thermal  transients.  Mode  buildup  transients 
occur  over  a  relatively  short  period  of  time  (10-60  ns)  that  depends  on  the 
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Copper  heat  sink 


Metallization 
Substrate  (n-GaAs) 
Confining  Layer  (GaAs) 
Graded  Cladding  (GaAIAs) 
Active  Layer  (SOW) 
Graded  Cladding  (GaAIAs) 
SiOj  current  isolation 
a  Cap  Layer  (p+GaAs) 
Metalization 

Solder  (In) 


0  448 

85.0 

0.5 

2.7 
0.01 

1.7 

0.24 

0.24 

3.0 


FIGURE  2.2  Broad-stripe  GRINSCH  SOW  laser  diode. 
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FIGURE  2.4  Output  waveforms  for  twelve  dominant  modes  of  the  SQW-60-300  laser 
modulated  at  105+100  mA  with  an  80  ns  pulse. 


laser  length.  When  the  laser  is  first  pulsed  from  threshold,  many  axial  modes 
are  excited  because  there  are  only  slight  differences  in  the  gain  factor  for  the 
different  modes  and  the  photon  population  in  each  of  the  modes  (which  drives 
the  stimulated  emission  rate)  is  due  initially  to  spontaneous  emission.^  But 
after  a  while,  the  slight  differences  in  the  gain  factor  cause  greater  and 
greater  differences  in  photon  populations  between  the  modes  because  a  larger 
photon  population  will  cause  a  larger  stimulated  emission  rate  for  that 
particular  mode.  In  other  words,  feedback  gradually  induces  most  of  the 
power  to  concentrate  in  a  few  central  modes  even  when  there  are  only  slight 
differences  in  their  gain  factors.  This  process  takes  approximately  10-60  ns 

depending  on  the  photon  lifetime  in  the  cavity  (which  is  related  to  the  laser 
length).  Thermal  transients  have  much  longer  time  constants  (several 
hundred  nanoseconds)  and  so  the  modal  instabilities  caused  from  these  can  last 
several  microseconds.  It  should  be  noted  that  the  severity  of  both  the  mode 
buildup  transients  and  the  thermal  transients  is  directly  related  to  the 
modulation  depth  of  the  pulse.  Figure  2.4  shows  that  for  the  SQW-60-300  laser, 
large  mode  buildup-time  transients  occur  over  the  first  10-20  ns  into  the  pulse 
and  then  the  pulse  stabilizes  out  to  a  fairly  flat  waveform.  Figure  2.5  shows  a 
laser  pulse  with  approximately  the  same  modulation  currents  as  for  Figure  2.4 
extended  out  to  1000  ns.  As  can  be  seen,  the  modes  maintain  fairly  flat 
waveforms  out  to  1000  ns,  indicating  negligible  thermal  transients.  A  rough 
picture  of  the  gain  spectrum  for  the  SQW-60-300  laser  can  be  derived  by  taking 
the  amplitude  of  each  of  the  modes  in  Figure  2.5  at  a  specific  point  in  time 
during  the  pulse  and  plotting  these  points  against  the  frequency  spacing 
between  the  modes.  (For  this  laser  the  modes  arc  spaced  approximately  0.3  nm 
apart  in  wavelength.)  It  can  be  seen  from  Figure  2.5  that  the  gain  spectrum 
for  the  SQW-60-300  laser  for  these  particular  operating  conditions  is  not  a 
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FIGURE  2.5  Output  waveforms  for  twelve  dominant  modes  of  the  SQW-60-300  laser 
modulated  at  108+100  mA  with  a  1000  ns  pulse. 


smooth  I.orenzian  shaped  curve  but  has  a  tew  ripples.  This  is  probably  due  to 
competition  between  modes  causing  some  of  the  modes  to  dominate  over  others. 
The  modes  of  the  80  ns  pulses  do  not  correspond  to  the  modes  of  the  1000  ns  and 
1600  ns  pulses  for  the  same  lasers  in  Figures  2.4  through  2.11.  This  is  because 
for  the  longer  pulses  the  temperature  in  the  active  region  at  the  start  of  the 
pulse  will  be  different  from  that  for  the  80  ns  pulse  because  for  the  longer 
pulse  the  temperature  has  more  time  to  drop  between  pulses. 

Figures  2.6  and  2.7  show  the  output  waveforms  for  the  itn  dominant 
modes  of  the  SQW-60-600  laser  modulated  at  a  50%  duty  cycle  irom  near 
threshold  to  100  mA  above  threshold  (and  producing  65  mW  of  light). 

Figure  2.6  shows  the  waveforms  up  to  80  ns  and  Figure  2.7  up  to  1000  ns.  For 
this  long  600  pm  length  laser,  it  is  Seen  that  the  mode  buildup-time  transients 
last  over  60  ns,  but  then  the  modes  stabilize  and  remain  fairly  flat  up  to 
1000  ns. 

Modui  instabilities  for  the  SQW-5-300  are  shown  in  Figures  2.8  and  2.9 
for  80  ns  and  1000  ns.  This  nairov^  stripe  laser  was  modulated  at  76-t-52  mA  at  a 
50%  duty  cycle,  producing  42  mW  of  light.  At  this  modulation  depth,  this  laser 
demonstrates  both  mode  buildup  transients  and  thermally-induced  transients. 
The  mode  buildup-time  transients  are  more  severe  and  last  10-20  ns,  whereas 
the  thermally-induced  transients  cause  gradual  changes  in  the  mode  powers 
which  last  over  1000  ns.  Figure  2.8  shows  that  the  shorter  wavelength  modes 
(beginning  ••’t  the  upper  left)  lose  power  during  the  pulse  and  the  longer 
wavelength  modes  gain  power  during  the  pulse.  From  this  figure  the  peak  of 
the  gain  curve  is  seen  to  shift  through  the  modes  as  time  progresses  and  the 
laser  heats  up  due  to  the  increased  current. 

The  ten-stripe  laser  (MQW- 100-250)  output  can  be  separated  into  several 
modes  for  each  of  the  ten  stripes.  Figure  2.10  shows  six  dominant  modes  shown 
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FIGURE  2.6  Output  waveforms  for  ten  dominant  modes  of  the  SQW-60-600  laser  modulated 
at  105+100  mA  with  an  80  ns  pulse. 
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FIGURE  2.7  Output  waveforms  for  ten  dominant  modes  of  the  SQW-60-600  laser  modulated 
at  101  +  100  mA  with  a  1000  ns  pulse. 
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FIGURE  2.9  Output  waveforms  for  six  dominant  modes  of  the  SQW-5-300  laser  modulated  at 
76+52  mA  with  a  1000  ns  pulse. 
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FIGURE  2.10  Output  waveforms  for  six  dominant  modes  of  the  10-stripe  laser 
diode  array  modulated  at  200  +  130  mA. 
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with  the  output  from  all  the  stripes  combined  for  each  mode.  The  pulse  length 
is  100  ns  with  a  50%  duty  cycle.  The  modulation  is  200  +  130  mA,  with  a  peak 
optical  power  of  90  mW. 

The  narrow-stripe  MQW  laser  (MQW-4-250)  modal  transients  are  shown 
over  a  periods  of  80  ns  and  1.6  ps  in  Figures  2.11  and  2.12.  The  modulation 
current  at  a  50%  duty  cycle  is  60+32  mA  for  a  peak  output  power  of  20  mW.  The 
thermal  transients  for  this  laser  are  particularly  severe  causing  modes  to 
completely  appear  and  disappear  over  a  period  of  a  few  hundred  nanoseconds. 

Temperature  changes  in  the  active  region  causes  the  positions  of  the 
modes  to  shift  as  well  as  the  relative  intensity  of  the  modes  to  shift.  These  two 
phenomena  are  respectively  known  as  mode  chirping  and  mode  hopping.  In 

the  case  of  mode  chirping,  the  actual  wavelength  of  the  individual  modes 
shifts.  In  the  case  of  mode  hopping,  the  gain  curve  shifts  causing  the  relative 
intensity  of  the  modes  to  change.  The  emission  wavelength  for  a  typical 
GaAlAs  laser  will  shift  an  average  of  3  This  is  partly  due  to  mode 

hopping  and  partly  due  to  mode  chirping.  Mode  hopping  is  the  dominant 
effect,  however  mode  chirping  was  measured  to  be  0.6  A/^C  by  Thompson^ 
and  0.93  A/°C  by  Dagenais,  et  al.,^  for  AlGaAs  lasers.  In  our  lab  we  measured 
2.5  A/°C  and  0.6  A/°C  for  the  total  modal  shift  and  for  mode  chirping  for  the 
SQW-60-300  laser.  For  the  lasers  with  large  thermally-induced  transients, 
mode  chirping  can  be  severe  enough  to  cause  the  mode  to  move  across  the 
detector  area  during  the  pulse  duration.  For  the  experimental  data  in  this 
report,  the  corresponding  detector  diameter  in  wavelength  units  was  0.4  A 
and  was  dependent  on  the  size  of  the  detector.  The  modes  had  roughly  a 
Gaussian  shape  with  a  FWHM  of  -1  A.  If  the  thermal  shift  is  greater  than 
-0.5  then  the  detector  will  shift  through  the  longitudinal  mode  during  the 
pulse,  and  therefore  the  data  will  not  accurately  measure  the  degree  of  mode 
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Tima, 


hopping.  For  the  MQW-4-250  and  the  SQW-5-300  lasers  at  the  long  pulse 
lengths  this  is  the  case.  However,  the  data  for  the  shorter  80  ns  excitation 
pulses  are  still  representative  of  the  degree  of  mode  hopping. 

One  way  to  measure  the  severity  of  the  modal  transients  in  a  10-stripe 
laser  is  to  measure  the  average  standard  deviation  of  the  ten  intensity 
waveforms  for  each  of  the  dominant  modes  from  a  square  wave  when  the  laser 
is  electrically  driven  by  a  square  current  pulse  above  the  bias  current.  This 
was  done  for  the  10-stripe  laser  at  a  constant  modulation  depth  of  130  mA  by 
varying  the  bias  current  around  the  threshold  value  of  200  mA.  Figure  2.13 
shows  the  results  plotted.  At  a  bias  of  190  mA,  slightly  below  threshold,  there 
is  a  standard  deviation  of  0.10  whereas  at  210  mA,  slightly  above  threshold, 
there  is  a  standard  deviation  of  0.17.  These  are  large  differences  in  the 
severity  of  the  modal  transients.  These  kind  of  fluctuations  are  also  apparent 
when  the  modulation  depth  is  changed.  In  general,  the  modal  transients  are 
less  severe  with  lower  modulation  depth. 
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FIGURE  2.13  Average  standard  deviation  from  a  square  wave  intensity 
profile  versus  the  bias  current  for  a  lO-stripe  laser  diode  array 
(modulation  depth  =  130  mA). 

These  results  can  be  qualitatively  explained  by  assuming  that  the 
position  of  the  gain  curve  with  respect  to  the  mode  positions  during  the  pulse 
plays  a  rrle  in  determining  the  severity  of  the  transients.  In  one  case  the 

gain  curve  is  centered  over  one  mode,  and  in  another  case  it  is  centered 
between  two  modes.  In  the  second  case  the  two  modes  may  compete  more 
strongly  for  power  resulting  in  a  longer  transient.  Gray  and  Roy^  have 
examined  mode  stability  in  AlGaAs  laser  diodes.  They  have  observed  these 
lasers  to  operate  in  different  regimes  of  single-mode  behavior(with  most  of 
the  power  concentrated  in  one  mode),  bistability  (where  two  modes  are  lasing 
but  only  one  at  a  time)  and  multimode  behavior  (with  both  modes  oscillating 
simultaneously)  depending  on  small  changes  in  operating  current  and 
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temperature.  They  utilize  a  multidimensional  laser  model  with  Langevin  noise 
sources  and  asymmetric  nonlinear  gain  to  model  this  behavior. 

2.2  General  discussion  of  experimental  results  obtained  using 
the  APD 

Our  experimental  data  show  that  two  types  of  modal  transients, 
mode-buildup  transients  and  thermal  transients,  occur  as  a  result  of  current 
injection  in  a  diode  laser.  The  mode-buildup  transients  occur  over  a  relatively 
short  period  of  time  (10-60  ns)  that  is  dependent  on  the  laser  length.  They  last 
10-30  ns  for  300  pm  length  lasers  and  20-60  ns  for  600  pm  length  lasers. 

Thermal  transients  are  distinguishable  from  mode  buildup  transients  by  much 
longer  time  constants  (on  the  order  of  several  hundred  nanoseconds)  and  so 
the  modal  instabilities  caused  from  these  can  last  several  microseconds. 

As  wc  have  seen,  the  broad-siripe  GRINSCH  SQW  laser  diodes  have 
negligible  thermal  transients  for  a  modulation  depth  of  100  mA.  The  MQW 
narrow  stripe  laser  had  the  largest  thermal  transients.  The  narrow-stripe  SQW 
laser  also  had  large  thermal  transients.  It  was  found  for  our  test  setup  that 
mode  chirping  can  cause  the  mode  to  move  out  of  the  detector  area  during  the 
pulse  duration  for  lasers  with  thermal  transients  greater  than  0.5  °C.  For  these 
lasers,  an  accurate  measure  of  the  mode  hopping  phenomena  can  only  be 
made  either  for  short  pulse  durations  or  at  low  levels  of  modulation. 

The  severity  of  both  the  mode-buildup  transients  and  thermal 
transients  is  directly  related  to  the  modulation  depth  of  the  pulse.  It  was  also 
found  that  the  severity  of  the  modal  transients  is  dependent  on  the  bias 
current  during  modulation. 


20 


2.3  Determination  of  the  relaxation  oscillation  frequencies  of  the 
laser  diodes  with  a  streak  camera 
The  relaxation  oscillation  frequencies  of  the  laser  diodes  as  a  function 
of  modulation  current  were  determined  with  a  streak  camera.  The  relaxation 
oscillations  in  the  laser  are  due  to  the  time  delay  involved  in  the  cyclical 
process  of  building  up  a  population  inversion  and  the  subsequent  stimulated 
emission  that  increases  the  photon  number  and  reduces  the  carrier  population 
to  below  the  inversion  level.  The  relaxation  oscillation  frequency  also 
determines  the  ultimate  modulation  frequency  at  which  a  semiconductor  laser 
can  be  modulated. 

Figure  2.14  shows  the  experimental  configuration  used. 


FIGURE  2.14  Streak  camera  experimental  configuration. 

With  this  setup,  the  laser  is  square  pulse  modulated  at  different  depths 
and  its  output  focused  onto  the  slit  of  the  streak  camera.  The  pulse  period  was 
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20  ns  with  a  duty  cycle  of  0.002%.  The  pulse  generator  had  a  rise  time  of 
150  ps.  In  a  streak  camera,  photons  enter  through  a  slit  and  are  focused  onto 
photocathode.  The  photocathode  converts  the  light  to  electrons  which  are 

accelerated  down  a  vacuum  tube,  optionally  amplified  with  a  microchannel 
plate,  and  then  deflected  across  a  phosphor  screen.  The  phosphor  screen  is 
read  out  with  a  CCD  camera.  The  streak  camera  has  the  ?dvantage  that  it  can 
observe  the  relaxation  oscillations  of  the  laser  diodes  with  picosecond 
resolution.  It  can  also  perform  both  single-pulse  and  multi-pulse 
measurements.  It  does  have  the  disadvantage,  though,  that  unless  the 
photocathode  obtains  fairly  large  amounts  of  light  the  streak  image  can  be 
noisy. 

Figures  2.15  and  2.16  show  the  streak  camera  plots  for  the  SQW-5-300 
laser  and  the  SQW-60-300  laser  taken  at  modulation  depths  of  80-t-143  mA  and 
0+254  mA,  respectively.  As  can  be  seen,  the  narrow  stripe  laser  gives  very 
clear  oscillations,  whereas  the  broad  stripe  laser’s  oscillations  are  hard  to  see 
and  quickly  become  chaotic.  The  broad  area  laser’s  tendancy  to  oscillate  with 
uncontrolled  filamentation  probably  contributes  to  this.  The  filamentation  is 
due  to  various  kinds  of  nonuniformity,  particularly  nonuniformities  in  the 
epitaxial  layers.^  This  type  of  behavior  has  been  studied  by  Yu,  et  al.,  in  the 
context  of  broad  area  operating  laser  arrays. Data  was  taken  for  all  four  of 

the  MacDac  lasers  for  a  certain  range  of  current  around  threshold.  The 
oscillation  frequency  versus  is  plotted  in  Figure  2.17.  As  predicted 

by  theory,^  ^  ^  ^  ^  a  straight  line  relationship  is  expected,  in  agreement  with 

our  experimental  results.  The  SQW-5-600,  SQW-5-600,  and  SQW-60-600  lasers 
were  pulsed  from  threshold  and  the  SQW-60-300  laser  was  pulsed  from  zero. 


FIGURE  2.15  Streak  camera  plot  of  SQW-5-300  laser  diode  modulated  at 
80  +  143  mA. 


Time,  ns 

FIGURE  2.16  Streak  camera  plot  of  SQW-60-300  laser  diode  modulated  at 
0  +  254  mA. 
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FIGURE  2.17  Oscillation  frequency  versus  normalized  current 
for  the  GRINSCH  SOW  laser  diodes. 
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These  data  can  be  compared  to  the  results  obtained  from  the  rate  equation 
modeling  in  Chapter  4. 
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CHAPTER  3 


THERMAL  MODELING  WITH  NASTRAN  FINITE  ELEMENT  CODE 

3.1  Setting  up  the  NASTRAN  finite  element  code 

The  NASTRAN  finite  element  code*^  is  a  software  program  for  modeling 
temperature  transients  in  different  materials  due  to  localized  heat  sources.  It 

was  adapted  to  simulate  temperature  transients  within  the  active  area  of 
semiconductor  lasers.  In  this  model,  the  diode  is  divided  into  a  grid  of  cells. 

The  conductivity,  density  and  heat  capacity  of  the  material  within  each  cell 
are  specified  as  well  as  any  heat  generated  from  photon  absorption, 
non-radiative  recombination,  or  ohmic  heating.  The  model  uses  the  thermal 

conductivities,  the  thermal  diffusivities  and  the  assumption  that  heat  can  flow 
out  only  through  a  constant  temperature  heat  sink  to  determine  the  thermal 
distribution  in  the  laser  during  modulation  or  in  steady-state  operation.^ 5, 16 
Both  2-dimensionai  and  3-dimensional  modeling  of  the  laser  diode  and 
heat  sink  were  used.  The  laser  itself  can  be  accurately  modeled 
two-dimensionally  because  edge  effects  are  assumed  not  to  be  important. 
Three-dimensional  effects  of  thermal  dissipation  in  the  heat  sink  become  a 

factor  as  the  heat  sink  thickness  is  increased  beyond  about  100  pm. 
Three-dimensional  modeling  is  very  computationally  intensive  so  it  was  used 

only  to  determine  the  thermal  impedance  of  the  heat  sink.  A  thickness  for  the 
two-dimensional  heat  sink  was  chosen  by  trial  and  error  until  the  temperature 
rise  in  the  heat  sink  matched  the  value  determined  three-dimensionally.  The 
transient  analysis  was  done  2-dimensionally. 

Three  cases  were  simulated  in  this  paper,  the  60  pm  wide  broad-stripe 
laser  (SQW-60-300),  the  5  pm  wide  narrow-stripe  laser  (SQW-5-300),  and  the 
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10-stripe  laser  diode  array  (MQW-100-250).  All  lasers  were  mounted 
junction -side-down. 

Figure  3.1  shows  the  2-dimensional  grid  used  in  the  simulation  for  the 
10-stripe  laser  array.  The  heat  sink  is  only  partially  shown  because  of  its  size. 
The  grid  was  built  with  variable  size  ceils  to  reduce  the  total  number  of  cells 
and  therefore  the  computational  time  required  to  run  the  simulation.  The 
first  55  horizontal  cells  are  each  1  pm  thick.  The  remaining  ceils  gradually 
increase  in  size.  The  active  area  of  the  laser  and  the  layers  with  large  heat 
fluxes  were  given  smaller  grid  sizes  so  the  temperature  transients  could  be 
more  accurately  modeled  in  those  regions.  Only  half  of  the  laser  was  modeled 
because  of  symmetry  about  the  vertical  axis.  The  loss  of  heat  into  the  air  was 
assumed  to  be  negligible.  The  ambient  temperature  was  fixed  at  the  bottom  of 
the  heat  sink  where  heat  was  allowed  to  flow  out  passively.  Figure  3.2  shows  a 
full  view  of  the  3-dimensional  model  used.  The  3-dimensional  model  used 
larger  cell  sizes  than  the  2-dimensional  model. 

Figure  2.3  shows  the  composition  and  layer  thicknesses  for  the  10-stripe 
laser  diode  (MQW-100-250).  The  metallization^^,  diode^ ,  and  heat 
sink 22  compositions  and  thicknesses  were  obtained  from  several  sources. 
Typically,  the  active  area  for  a  MQW  laser  consists  of  four  0.013  pm  thick  layers 
of  Ga.94Al.06As  separated  by  three  .004  pm  thick  layers  of  Ga.8A1.2As'^.  In  this 
work  the  active  area  was  modeled  as  a  single  layer  0.06  pm  thick  with  the 
average  composition  of  Ga.91  AI.09AS.  The  ten  stripes  were  defined  by  proton 
implantation.  The  heat  capacity  and  conductivity  of  GaAs  were  used  for  the 
proton-implanted  pGaAs  (p+GaAs)  current-isolating  layer.  The  heat  capacities 
of  Ga.oAl  4AS  and  Ga.91Al.09As  were  not  available  and  were  approximated  by 
doing  a  linear  fit  of  data  between  the  heat  capacity  values  of  GaAs  and 
Ga.5Al.5As.  The  2-dimensionaI  heat  sink  thickness  used  was  300  pm,  which 
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FIGURE  3.1  Two-dimensional  finite  element  grid  structure  for  thermal 
modeling  of  laser  diode  array 


FIGURE  3,2  Three-dimensional  grid  structure  for  thermal  modeling  of 
laser  diode  heat  sink. 
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approximately  matched  the  thermal  impedance  of  the  3-dimensional  heat  sink. 
The  width  of  the  heat  sink  was  600  pm.  The  temperature  of  the  heat  sink  on 
the  bottom  row  of  cells  was  fixed  at  20®C. 


Figure  2.2  shows  the  composition  and  layer  thicknesses  for  the 
broad-stripe  laser  diode  (SQW -60-300).  The  narrow  stripe  SQW-5-300  laser  diode 
has  the  same  structure  but  with  a  S  pm  stripe  width  instead  of  a  60  pm  stripe 
width.  The  stripes  are  designated  by  a  Si02  layer.  The  thermodynamic 
properties  of  the  materials  present  in  the  lasers  are  given  in  Table  3.1. 

TABLE  3.1 

Thermodynamic  Properties 


Material 

Molecular 

Weight 

Density 

Conductivity 

Heat 

Capacity 

Heat 

Capacity 

Thermal 

Diffusivity 

g/gmole 

g/cm^ 

W/(pm-OC) 

J/(OC-pm3) 

cm2/sec 

Ni 

58.71 

8.9023 

O.9IE-423 

6.2223 

395.E-14 

In 

114.82 

7.3023 

O.8I8E-423 

6.4323 

171.E-14 

Au 

196.97 

18.8823 

EHISSH 

6.0723 

243.E-14 

Cr 

52.00 

7.2023 

0.939E-423 

5.5623 

322.E-14 

Sn 

118.69 

6.5423 

0.617E-424 

6.4523 

149.E-14 

Cu 

63.54 

8.9223 

4.OIE-422 

5.8423 

343.E-14 

Ti 

47.90 

4.523 

0.219E-423 

5.9823 

235.E-14 

Ge 

72.59 

5.3523 

0.602E-423 

5.5823 

172.E-14 

Pt 

195.09 

21.4523 

O.7I6E-423 

6.1823 

284.E-14 

W 

183.85 

19.3523 

1.73E-423 

5.8023 

255.E-14 

GaAs 

144.64 

5. 3125 

0.45E-426 

11.0523 

170.E-14 

0.2727 

Ga.5Al.5As 

- 

0.12E-426 

- 

235.E-14 

0.05128 

- 

0.12E-426 

- 

222.E-14 

Ga.4Ai.6As 

- 

0.12E-426 

- 

248.E-14 

Ga.7Ai.3As 

- 

0.13E-426 

- 

209.E-14 

Ga.85Ai.15As 

- 

O.I6E-426 

- 

190.E-14 

Ga.9i  A1  09AS 

- 

O.2IE-426 

- 

182.E-14 

Ga 

69.72 

5.90423 

- 

6.20523 

220.E-14 

A1 

26.98 

2.70223 

2.I8E-424 

5.8223 

244.E-14 

As 

74.922 

5.7224 

0.502E-423 

5.8823 

188.E-14 

Si02 

60.08 

2.3223 

O.OlE-422 

10.824 

174.E-14 

TO 
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Determination  of  heat  sources  within  the  diode 


Heat  is  generated  in  the  laser  diode  from  ohmic  heating  at  the 
p-contact,  non-radiative  recombination  in  the  active  layer,  photon  absorption 
in  the  active  layer,  and  photon  absorption  in  the  substrate  and  p-cap 
layer. 22, 29, 30,6, 12 

The  equations  in  Table  3.2,  derived  in  a  paper  by  Kobayashi  and 
Furukawa^®,  are  used  to  calculate  the  heat  input  to  the  various  layers  of  the 
diode. 

TABLE  3.2 

Equations  for  Calculating  Heat  Input  in  the  Laser  Diode 


below  threshold 

(Kith) 

above  threshold  (l>Ith) 

Layer 

VjI(l*'nspon) 

+  Vjl  TisponC  1  -0 

( Vjl- VjIih)(Tii-Tiext)''''^jIthhspon(I  *f) 
+Vjlth(l  'nspon)+(  Vjl-VjIth)(l-Tispon)'n’ 
+(VjI-V  jIth)(l*Oh'hsDon 

active  layer 

V  jlthh  st)on(f/2)+(ViI- V  iIth)'n''nsoon(f/2) 

GaAs  substrate 

V  iIthhsDon(f/2)+(V  jl-V  iIth)'n'hsDon(f/2) 

p-GaAs  cap  layer 

i2Rs 

I^Rs 

each  layer  except 
active  layer 

In  these  equations,  Vj  is  the  junction  voltage,  Tl  i  is  the  internal 
efficiency  of  lasing,  Rs  is  the  series  resistance  of  each  layer,  T|spon  is  the 

internal  quantum  efficiency  of  the  spontaneous  emission,  Tl’  =  1-  Tli,  Ith  is  the 
threshold  current,  T|exi  is  the  external  differential  quantum  efficiency  of 

lasing,  I  is  the  current,  and  f  is  the  fraction  of  spontaneous  emission  radiated 

from  the  active  layer  and  absorbed  by  the  GaAs  substrate  and  p-GaAs  cap  layer. 

The  value  f  is  approximately  equal  to  1  -  (l-(ni/no)^)-^  ,  where  Uq  and  nj  arc 
the  refractive  indices  of  the  active  layer  and  the  cladding  rc.spcctivcly.2  2 
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For  T|i=l  (T1'=0)  and  Tlspon=l>  assuming  ohmic  heating  to  occur  primarily 
at  the  p-contact  metallization  layer^^^  the  equations  reduce  to  those  shown  in 
Table  3.3. 


TABLE  3.3 


Simplified  Equations  for  Calculating  Heat  Input  in  the  Laser  Diode 


at  threshold  (I=Ith) 

above  threshold  (l>Ith) 

Layer 

Viith  (1-0 

(ViI-ViIth)(l-i}ext)+Vilth(l-f) 

active  layer 

Vi  Ith  f/2 

ViIth(f/2) 

GaAs  substrate 

Vi  Ith  f/2 

ViI,h(f/2) 

p-GaAs  cap  layer 

Ith^Rs 

i2Rs 

p-contact  layer 

The  operating  parameters  for  the  diodes  are  given  in  Table  3.4.  The 
operating  parameters  for  the  MQW  lasers  were  obtained  from,  their  data  sheets. 
The  operating  parameters  for  GRINSCH  SQW  lasers  were  obtained  from 
D.K.  Wagner,  et.  al.^  and  data  sheets.  V{,j  is  the  diode  voltage  at  threshold,  and 
Vb2  is  the  diode  voltage  at  the  operating  current  above  threshold.  The 
narrow-stripe  laser  has  a  higher  ohmic  resistance  because  of  the  smaller 
contact  area. 


TABLE  3.4 


Operating  Parameters  for  Laser  Diodes 


Parameter 

SQW-60-300 

MQW-100-250  |SQW-5-300 

Laser  type 

broad-stripe 

10-stripe  1-stripe 

Optical  power  (mW) 

80 

200 

100 

200 

46 

Vi(V) 

1.51 

1.51 

1.51 

1.51 

1.51 

Vbl  (V) 

1.54 

1.54 

1.76 

1.76 

- 

Vb2  (V) 

1.63 

1.78 

1.81 

1.9 

- 

Ith  (mA) 

100 

100 

200 

200 

75 

1  (mA) 

200 

340 

340 

475 

128 

Rs  (G) 

.6 

.6 

.85 

.85 

6 

f 

.65 

.65 

.65 

.65 

.65 

hext 

.473 

.473 

.49 

.49 

.62 
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When  the  operating  parameters  of  Table  3.4  are  substituted  into  the 
equations  of  Table  3.3,  the  generated  heat  in  the  specified  layers  of  the  laser 
can  be  found.  These  values  are  listed  in  Table  3.5.  The  listed  values  for  total 
heat  generated  represent  the  total  power  into  the  laser  minus  the  optical 
power  output.  For  a  10-stripe  laser,  each  stripe  will  receive  only  one  tenth  of 
these  power  levels.  The  broad-stripe  lasers  are  more  efficient  than  the 
10-stripe  array  because  current  spreading  is  not  as  significant  for  them.  It 

occurs  only  at  the  two  outer  edges  of  the  broad  active  area  as  opposed  to  on 
either  side  of  each  stripe  for  the  10-stripe  laser. 

TABLE  3.5 

Heat  Generated  in  the  Laser  Diode  (in  Watts) 


Parameter 

SQW-60-300 

MQW-100-250  1  SQW-5-300 

Laser  type 

broad-stripe 

10-stripe  1-stripe 

Optical  power  (mW) 

80 

200 

100 

200 

46 

at  threshold: 

active  layer 

0.0529 

0.0529 

0.106 

0.106 

0.040 

GaAs  substrate 

0.049 

0.049 

0.098 

0.098 

0.037 

p-GaAs  cap  layer 

0.049 

0.049 

0.098 

0.098 

0.037 

p-contact  layer 

0.006 

0.006 

0.034 

0.034 

0.034 

above  threshold; 

active  layer 

0.132 

0.244 

0.214 

0.317 

0.070 

GaAs  substrate 

0.049 

0.049 

0.098 

0.098 

0.037 

p-GaAs  cap  layer 

0.049 

0.049 

0.098 

0.098 

0.037 

p-contact  layer 

0.024 

0.069 

0.098 

0.192 

0.098 

Total  heat  generated 

at  threshold: 

0.157 

0.157 

0.336 

0.336 

0.147 

above  threshold: 

0.255 

0.41  1 

0.508 

0.706 

0.241 

For  the  MQW-100-250  and  SQW-5-300  lasers  the  current  in  the  active 
layer  spreads  out  laterally,  increasing  the  effective  stripe  width. 

Determination  of  the  actual  current  distribution  in  the  diode's  active  area  was 
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accomplished  by  powering  the  laser  below  the  threshold  current  and  imaging 
the  laser  facet  onto  a  CCD  camera.^^  The  CCD  camera  is  interfaced  with  an 
image  processor  which  can  plot  the  intensity  distribution  of  the  spontaneous 
emission.  The  assumption  was  made  that  the  output  intensity  distribution  is 
proportional  to  the  current  distribution.  For  the  MQW- 100-250  array,  with  cell 
widths  of  1  pm  the  fraction  of  power  per  cell  was  measured  to  be  0,  0.03,  0.10, 
0.17,  0.20,  0.20,  0.17,  0.10,  0.03,  0  across  the  ten  active  layer  cells  for  a  current 
10  mA  below  threshold.  This  profile  was  also  used  as  an  approximation  for  the 

substrate  and  cap  layers  For  the  broad-stripe  laser,  because  the  stripe  is  wide 

and  the  central  r  .i  of  the  stripe  was  the  area  of  interest,  no  current 
spreading  was  '  ..wd  in  the  simulation.  For  the  SQW-5-300  narrow-stripe  laser  a 
Gaussian  distribution  was  seen  with  a  full  width  half  maximum  (FWHM)  of 
20  M‘«  at  a  current  13%  of  threshold.  This  FWHM  decreased  as  the  current 

increased,  with  a  FWHM  of  15  pm  for  a  current  31%  of  threshold,  a  FWHM  of  10 

pm  for  a  current  44%  of  threshold,  a  FWHM  of  7  pm  for  a  current  67%  of 

threshold,  and  a  FWHM  of  5.5  pm  for  a  current  91%  of  threshold.  The  current 

distribution  above  threshold  cannot  be  accurately  measured  by  imaging 
because  the  lasing  regions  will  show  disproportionately  more  intensity  than 
the  non-lasing  regions  because  the  lasing  light  is  directional.  For  this  work, 
the  value  of  5  pm  for  the  FWHM  was  used  to  model  the  SQW-5-300  narrow-stripe 
laser. 
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3.3  Steady-state  results 

Table  3.6  shows  the  change  in  temperature  of  the  active  area  in  the 
center  of  the  stripe  from  the  heat  sink  temperature  of  20  °C  for  the  three  laser 
diodes  for  various  current  levels  run  at  steady- state. 


TABLE  3.6 

Change  in  Temperature  of  Active  Area  during  Steady-State  Operation 


Laser 

Temperature  Rise,  oC 

SQW-60-300 

80 

4.78 

200 

7.60 

MOW-100-250 

200 

10.7 

SOW-5-300 

46 

10.24 

3.4  Transient  results  and  discussion 

Table  3.7  lists  the  temperature  change  of  the  active  area  in  the  center 
of  the  stripe  for  the  three  lasers  during  modulation  with  80  ns  and  1600  ns 
pulses.  For  the  80  ns  pulse,  a  50%  duty  cycle  was  used.  For  the  1600  ns  pulse, 
the  laser  was  set  initially  at  the  threshold  current  level  until  the  temperature 
reached  steady  state  and  then  the  modulating  current  was  applied.  This 
approximates  a  50%  duty  cycle  pulse  since  in  1600  ns  the  temperature  in  the 
active  region  will  have  a  long  enough  time  to  settle  back  to  the  threshold 
temperature  after  the  modulated  pulse  ends.  The  transient  simulations  were 
run  with  varying  step  sizes.  The  first  section  was  run  with  a  step  size  of  9  ps  to 
bring  the  diode  near  steady  state.  After  that,  a  step  size  of  1  ns  was  used. 
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TABLE  3.7 


Thermal  Transients  in  Active  Area  of  Pulsed  Laser  Diodes 


Laser 

Peak  Power,  mW 

Temperature  Rise,  C 
(80  ns  pulse) 

Temperature  Rise,  C 
(1600  ns  pulse) 

SQW-60-300 

80 

.092 

.58 

200 

.219 

- 

MQW-100-250 

200 

0.35 

- 

100 

0.18 

- 

SOW-5-300 

46 

0.396 

2.12 

The  SQW-60-300  broad-stripe  laser  has  smaller  transients  than  the 
MQW- 100-250  ten-stripe  array  because  it  has  minimal  current  spreading.  It 
also  has  a  larger  total  active  area  (60  pm  versus  40  pm  for  the  10-stripe  laser), 
and  has  a  longer  length  of  300  pm  versus  -250  pm  for  the  10-stripe  laser. 
Figure  3.3  plots  the  temperature  in  the  active  area  of  the  laser  during 
modulation  for  the  SQW-60-300  laser  for  a  peak  optical  power  of  80  mW  and 
80  ns  pulses.  Figure  3.4  shows  the  SQW-60-300  laser  with  a  peak  optical  power 
of  80  mW  for  a  1600  ns  pulse. 

These  temperature  transients  can  be  translated  into  modal  transients 
and  compared  to  experimental  data  by  using  the  relationship  that  the  emission 
wavelength  for  a  typical  GaAlAs  laser  will  shift  an  average  of  -3  A/°C.‘^ 

Table  3-8  shows  the  thermal  time  constants  of  the  transients  fitted  to  the 
equation 

AT(t)  =  (kl(l-e(-‘/Tl))+k2(]-e(-‘/T2)))  (3.1) 

where  AT(t)  is  the  temperature  rise,  t  is  the  time  in  ns,  and  kl  and  k2  are 
steady  state  temperature  values  associated  with  the  time  constants  T1  and  T2. 
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Time,  ns 

FIGURE  3.3  Thermal  transients  in  active  area  of  SQW-60-300  laser  during 
modulation  with  a  160  ns  square  wave  and  80  mW  peak  optical 
power. 
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FIGURE  3.4  Thermal  transients  in  active  area  of  SQW-60-300  laser  during 
modulation  with  a  1600  ns  pulse  -ad  80  mW  peak  optical  power. 


TABLE  3-8 

Thermal  Time  Constants  Fitted  to  Transient  Data  from  the  Finite  Element  Code 

Analysis 


Laser 

Transient 
Time,  ns 

Output 
Power,  mW 

n 

IHH 

SQW-60-300 

80 

0.0281 

4.5 

0.091 

■BB 

1600 

0.295 

400 

0.41 

■BB 

SQW-5-300 

80 

0.16 

■BB 

0.45 

110 

FWHM=5 

1600 

0.8 

2.1 

1600 

3.5  Effect  of  varying  cell  sizes 

The  cell  size  was  varied  for  the  lasers  to  determine  its  effect  on  the 
simulations.  In  one  case,  cell  widths  of  1  |j.m  were  used  with  cell  thicknesses  of 
0.5  tim  in  the  confinement  area.  In  the  other  case  cell  widths  of  0.1  jim  were 
used  with  cell  thicknesses  of  0.05  in  the  confinement  area.  Small  changes 
in  the  thermal  transients  on  the  order  of  a  few  percent  occurred,  indicating 
that  accurate  simulations  can  be  made  with  the  larger  cell  sizes. 
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CHAPTER  4 

TEMPERATURE-DEPENDENT  RATE  EQUATIONS  MODEL 


4.1  Introduction 

Rate  equations  model  the  modal  output  power  characteristics  of  laser 
diodes.  The  characteristics  and  dimensions  of  the  laser  are  entered  into  the 
program  along  with  the  driving  current,  and  the  equations  are  solved 
numerically  to  predict  the  power  in  each  of  the  laser  modes  as  a  function  of 
time.  The  rate  equations  model  is  used  to  obtain  information  on  relaxation 
oscillation  frequency,  modal  power  distribution,  total  output  power,  and  the 
settling  times  of  the  modal  power  fluctuations  due  to  mode-buildup  time 
transients  and  thermal-induced  transients. 

Temperature-independent  rate  equations  have  been  described  by  many 
workers  in  the  field  including  Marcuse  and  Lee^^,  Yariv32,  Adams  and 
Osinski33,  and  Agrawal  and  Dutta^^.  Byrne  and  Keating  have  developed 
temperature-dependent  rate  equations  by  adding  temperature  dependence  to 
the  threshold  current  density  and  to  the  peak  of  the  modal  gain  curve 

distribution,  however,  they  only  investigated  dc  behavior.35  xhis  report  uses 

the  equations  developed  by  Byrne  and  Keating  and  extends  their  work  by 
examining  transient  behavior. 

4.2  Description  of  Rate  Equations 

The  electron  and  photon  transient  behavior  in  laser  diodes  can  be 
described  by  the  following  rate  equations. 
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(4.1) 


dn  e _  i  He 
dt  ~  ed  Te 


n/2 

M  A  DD(ne- no)S'u 
x>=-n/2 


,  —  Dy  lie"  Su  +  MrA  DuCng-  no)S'u  (4.2) 

ai  Tr  tp 

The  variables  in  these  equations  are  defined  as  follows: 

ne  =  injected  electron  density  in  the  conduction  band, 

t  =  time  coordinate,  s 

J  =  current  density,  A/cm^ 

e  =  electron  charge  (1.6  x  10'^^  C) 

d  =  thickness  of  active  laser  region,  tim 

Xe  =  carrier  lifetime,  s 

Xr  =  carrier  lifetime  due  to  radiative  recombination  processes 

n  +  1  =  number  of  modes  with  gain 

A  =  stimulated  emission  factor,  pm^/s 

I>o  =  line  shape  factor 

no  =  carrier  density  for  transparency,  pm'^ 

So  =  photon  density  of  \)th  laser  mode,  pm"^ 

r  =  ratio  of  active  layer  volume  to  optical  mode  volume  (mode 

confinement  factor) 

y  =  fraction  of  spontaneous  emission  coupled  into  the  lasing 

mode  (assumed  to  be  frequency  independent) 
(spontaneous  emission  factor  ) 

Xp  =  photon  lifetime,  s 

M  =  gain  saturation  term 
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The  photon  lifetime,  Tp,  is  related  to  the  group  refractive  index  of  the 
laser  medium,  ng,  and  the  effective  cavity  loss  coefficient,  a,  by  the  relation 

Tp  =  ng/(c  a)  .  (4.3) 


The  photon  lifetime  represents  the  average  time  the  photon  spends  in 
the  cavity.  It  takes  into  account  cavity  losses  due  to  absorption  of  photons  in 
the  cavity  and  mirror  losses.  The  effective  cavity  loss  coefficient,  a,  is  related 
to  the  cavity  length  and  reflectivity  of  the  mirrors,  R,  by 


1  I  ^ 

a  =  tti  +  -  In  ^  , 


(4.4) 


where  R  =  V  R1*R2  and  R1  and  R2  are  the  intensity  reflectivity  coefficients 
(reflectances)  for  the  front  and  back  facets.  The  reflectance,  Rl,  is 
approximately  0.32  for  an  uncoated  facet. ^  The  intrinsic  loss  coefficient  is  a,. 
The  gain  compression  term,  M,  is  defined  as 


M  = - 

1  + 


_1 _ 

n/2 

v=-n/2 


(4.5) 


where  £  is  the  gain  compression  factor.  The  main  effect  of  the  gain 
compression  term  is  to  dampen  the  relaxation  oscillations  of  the  laser.  The 

gain  compression  term  (sometimes  referred  to  as  gain  saturation  factor)  is  a 
phenomenological  term  which  takes  into  account  various  nonlinearities 
occurring  in  the  laser  medium  and  spectral  hole  burning. It  has  also 
been  used  to  represent  the  effect  of  lateral  carrier  diffusion  in  narrow  width 
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lasers. The  transient  response  is  not  very  sensitive  to  the  exact  value  of 
the  gain  saturation  factor  (other  than  damping  the  relaxation  oscillations). 
Values  for  £  could  not  be  derived  accurately  from  the  experimental  data  in  this 
report,  though,  in  Chapter  5,  approximate  values  for  £  are  found  which 
moderately  dampen  the  relaxation  oscillations  for  the  various  lasers. 


Equation  (4.1)  can  be  normalized  by  multiplying  through  each  side  by 
TpA  as  shown  in  equation  (4.6).  Equation  (4.2)  is  normalized  by  as  shown 
in  equation  (4.7).  The  assumption  is  made  that  is  approximately  equal  to  Xe 
as  will  be  discussed  later  in  Section  4.10. 


d(neXpA)  jXpA  (neXpA) 

dt  ~  ed  Te 


^  A  ivusuXeA; 

X  D^((neXpA)-(tioXpA)}  - (4.6) 

v=-n/2  ^ 


d(Sa)XeA)  py 
dt  “  Xp 


D'u(netpA) 


xp  EDy { (neXpA)-(noXpA)} 


M(S'^XeA ) 


(4.7) 


Using  the  equalities 
N  =  ngXpA, 

Su  =  SyXgA, 

No  =  noXpA, 

_  JtcXpA 

""  ed  ’ 

we  obtain  the  following  equations. 


dt 


J  n 


n/2  \ 

N  -  M  Do(N-No)Su 
v=-n/2  j 


(4.8) 
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(4.9) 


dSv)  1  ^ 

—  =  ;^(rYDu  N  -  Sy  +  M  r  Dy(N-No)Sy). 

The  normalization  of  the  electron  and  photon  density  equations  not 
only  serves  to  simplify  the  equations  but  also  facilitates  the  use  of  a 
Runge-Kutta  numerical  analysis.  This  is  because  the  stiffness  of  the  equations 

is  determined  by  the  ratio  of  the  independent  variables.  The  ratio  of  the 
electron  density  over  the  photon  density  is  normally  quite  high  so  the 
equations  are  stiff  and  require  a  very  small  time  step  to  converge.  The 
normalization  factors  reduce  the  ratio  of  the  electron  and  photon  density 
terms  and  so  the  system  converges  with  a  much  larger  time  step. 

The  first  term  in  (4.1)  is  the  increase  in  electron  carriers  due  to 
current.  The  second  term  is  the  decrease  in  electron  carriers  due  to 
spontaneous  emission,  and  the  third  term  is  the  decrease  or  increase  in 
electron  carriers  due  to  stimulated  emission  or  absorption  of  photons. 

Equation  (4.2)  represents  the  change  in  photon  population  for  each 
separate  guided  mode.  The  first  term  in  (4.2)  is  the  increase  in  the  photon 
population  due  to  spontaneous  emission.  The  factor  y  accounts  for  the  fact  that 
only  a  fraction  of  the  spontaneously  emitted  photons  contribute  to  the  guided 
mode.  The  second  term  represents  the  loss  of  photons  due  to  cavity  losses.  The 
third  term  represents  the  stimulated  emission  or  absorption  of  photons.  Note 
that  the  stimulated  emission  rate  represented  by  this  term  is  due  both  to  a  gain 
factor  and  to  the  concentration  difference  of  electrons  above  transparency 
and  the  concentration  of  photons  in  that  particular  mode.  As  described  in 
Section  2.1,  it  is  the  fact  that  the  stimulated  emission  rate  is  driven  by  the 
concentration  of  photons  in  a  particular  mode  that  causes  the  mode-buildup 
transients.  It  is  only  because  of  feedback  that  the  power  tends  to  gradually 
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concentrate  in  a  few  dominant  modes.  The  first  and  third  terms  of  (4.2)  are 
multiplied  by  F  to  account  for  the  different  volumes  occupied  by  the  electrons 
and  photons.  This  can  be  verified  by  holding  F  constant  and  iterating  the  rate 
equations  for  different  values  of  active  layer  thicknesses.  The  optical  power 
out  of  the  cavity  is  independent  of  the  active  layer  thickness  as  long  as  F  stays 
constant. 

The  photon  density  can  be  converted  into  optical  power  out  of  the  laser 
cavity  as  shown  in  the  following  equation: 

n /2 

P-  X  ^-7^^  n„/(c('a-a.))  <"■'<» 


where 


E  = 


h  c 

Xpo 


c  id  the  variables  are  defined  as 
P  =  optical  power,  W 

S  =  width  of  laser  stripe,  pm 

L  =  cavity  length,  pm 

E  =  energy  per  photon,  J 

h  =  Planck’s  constant  (6.626  x  10"^’*  J/s) 
c  =  velocity  of  light  in  vacuum  (3  x  lO'*^  pm/s) 

^po  =  wavelength  at  the  peak  of  the  mode  distribution  at  time  0 

o-aj  =  facet  loss  =  ^  R  ‘ 

In  (4.10),  the  photon  density  is  multiplied  by  the  mode  volume  and 
energy  per  photon  to  give  the  energy  in  joules  represented  by  all  the  photons 
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in  the  cavity.  This  energy  is  then  divided  by  the  photon  lifetime  (without 
internal  losses)  to  give  the  optical  power  coming  out  of  the  cavity. 

4.3  Determination  of  the  stimulated  emission  factor  and  the 

transparency  electron  current  density 

The  stimulated  emission  factor.  A,  affects  the  oscillation  frequency  of 
the  relaxation  oscillations.  A  larger  stimulated  emission  factor  increases  the 
oscillation  frequency.  The  stimulated  emission  factor  can  be  obtained  through 
values  obtained  for  the  gain.  In  Thompson^  and  Kressel  and  Butler^  2  the 
theoretical  relation  between  peak  optical  gain,  g,  in  cm"',  and  injected 
current  density,  J,  in  Amps  cm*^,  for  pure  GaAs  at  various  temperatures  is 

g=P(J-Jo)  .  (4.11) 


where  Jo  is  the  transparency  electron  current  density.  The  gain,  g,  is  defined 
as  the  power  emitted  per  unit  volume  divided  by  the  power  crossing  per  unit 
area.  If  this  relation  is  plotted,  the  slope  of  the  line  represents  the  gradient 

P  (sometimes  referred  to  as  the  differential  gain  coefficient)  in  cm  Amp"'.  The 
relationship  between  P  and  the  stimulated  emission  factor.  A,  is  given  by 


^  _  P  d  c  c 

Te  ng 


(4.12) 


The  transparency  electron  density  for  the  peak  mode  of  the  distribution 
is  related  to  the  transparency  electron  current,  Jo,  by  the  relation ^  2 


A 

no 


Jo 

ed 


(4.13) 
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The  actual  threshold  current  of  the  device  can  be  determined  by  the 


threshold  condition  where  the  gain  equals  the  loss^. 


ln(l/R) 

rg=  Oii+—\ — \ 


(4.14) 


and  the  gain  current  relation  (4.11).  The  expression  derived  from  these  two 
equations  is 


((Jo 


a; 

r 


L  + 


ln(l/R) 

r  p 


(4.15) 


where  1th  is  the  threshold  current  during  the  pulse  in  Amps. 

In  narrow  stripe  gain-guided  laser  diodes  the  current  will  not  remain 
confined  to  the  narrow  stripe  width  defined  by  the  oxide  or  proton 
implantation.  There  are  two  processes  by  which  this  occurs.  First,  the  current 
will  spread  out  during  its  trip  from  the  metal  contact  down  to  the  active  area. 
This  effect  is  directly  related  to  the  current  density,  and  the  thickness  and 
conductivity  of  the  confining  layer  which  comes  between  the  active  layer  and 
the  metallization.  Once  the  current  reaches  the  active  layer,  it  will  diffuse 
outward  along  the  active  region.  This  effect  is  related  to  the  diffusivity  and  to 
the  carrier  lifetime  in  the  active  region. 

Casey  and  Panish40  have  derived  equations  that  predict  the  carrier 
distribution  in  the  laser  diode  active  layer  based  on  both  current  spreading 
and  active-layer  carrier  diffusion.  They  combine  the  current  spreading 
expressions  described  by  Yonezu  et.  al-^^  with  the  carrier-diffusion  model  of 
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Hakki'^2.  Thompson^  and  Streifer,  et.al.,'*^  also  describe  the  effects  of  current 
spreading  and  diffusion. 

In  the  Casey  and  Panish  approach,  current  spreading  effects  are 
characterized  by  a  current  spreading  length, 


‘0  - 


2  L 
K  Ps  lo 


(4.16) 


where 


(l+(SKps/(2L))lt)-^  -1 
Skps/(2L) 


1  1  _ L _ 

Ps~(Pl/dl)  (p2/d2) 


Here  Iq  is  the  lateral  spreading  current.  It  is  the  total  current,  Ps  is  the 
composite  sheet  resistivity,  di  and  d2  are  the  lengths  of  the  layers  through 
which  the  current  flows  to  reach  the  active  layer  with  their  respective 
resistivities,  pi  andpi,  and  k  is  the  exponential  junction  parameter  e/nkT.  In 
the  SQW  lasers  used  in  this  report,  the  confining  layer  on  the  p-side  is 
Ga.6Al.4As.  Using  a  mobility,  pp,  of  250cm2/(V  s)  44_  and  hole  doping,  p,  of 
2x10'^  cm'3,  this  give  a  conductivity  (g=Pp  p  e)  of  8  mho  cm'^.  The 
resistivity  is  the  inverse  of  the  conductivity.  For  a  thickness,  di  +  d2.  of  1.7  pm, 
stripe  width,  S,  of  5  pm  and  stripe  length,  L,  of  300  pm,  this  makes  a 
characteristic  spreading  length  of  1.2  pm  for  a  current  density  of 
8500  A/cm^,  1.64  pm  for  a  current  density  of  5000  A/cm^  (threshold  for  the 
SQW-5-300  narrow  stripe  laser),  2.12  pm  for  a  current  density  of  3330  A/cm^, 
and  6.68  pm  for  a  current  density  of  670  A/cm^.  Also,  it  should  be  noted  that 
for  proton-bombarded  stripe-geometry  lasers  (such  as  the  MQW-4-250  laser) 
lateral  current  spreading  is  further  decreased.'^  ^ 
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Diffusion  effects  in  the  active  layer  are  characterized  by  a  diffusion 
length.  The  diffusion  length  can  be  calculated  knowing  the  electron  mobility, 
He,  and  the  electron  lifetime  in  the  p-type  material,  Xp: 


where 


De^ 


PekT 

e 


(4.17) 


is  the  diffusion  constant,  k  is  the  Boltzman  constant,  1.38  10'23  j/k,  and  T  is  the 
temperature  in  degrees  Kelvin. 46  Por  electrons  in  GaAs  at  concentrations 
below  10^^  cm'^,  Pe  *s  0.85  m^/(V  sec). 46. 12  (This  value  will  decrease  at 
higher  electron  concentrations.)  Using  this  value,  a  temperature  of  300  °K, 
and  an  electron  lifetime  of  2  ns'  * ,  a  value  of  6.6  pm  is  calculated  for  the 
diffusion  length.  Hakki42  and  Streiffer43  use  values  of  5  and  3  pm , 
respectively,  for  the  diffusion  length  when  modeling  current  diffusion  in 
GaAs  lasers. 

For  stripe  widths  under  20  microns,  electron  spreading  and  diffusion 
effects  start  playing  a  major  role  in  the  mode  size  and  efficiency  and  so  values 
of  P  and  Jo  determined  from  wide  width  laser  diodes  are  no  longer  valid.  The 
effect  of  current  spreading  and  diffusion  on  the  gain  factors  can  be 
determined  by  first  predicting  the  profile  of  the  current  along  the  stripe. 

Then  new  effective  values  of  P  and  Jo  for  the  case  of  narrow  stripe  lasers  can 
be  calculated.  Casey  and  Panish  predict  the  electron  concentration  in  the 
active  layer.  Their  equations  can  be  normalized  to  give  the  following  two 
equations; 
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'  ('  -i^)  '  ^°'  '*' '  2  • 


Jr(x) 


for  Ixl  >  ~ 


where  Jr(x)  is  the  relative  percentage  of  current  density  at  any 
from  the  center  of  the  stripe. 

Figure  4.1  shows  Jr(x)  versus  x  for  a  stripe  width,  S,  of  5  pm, 
length,  Ln,  of  5  microns  and  a  spreading  length,  Iq,  of  1.2  pm. 


FIGURE  4.1  Relative  current  density  versus  distance. 


(4.18) 

(x-S/2)>> 
Ln  ) 

(4.19) 

distance  x 

diffusion 
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Equations  (4.18)  and  (4.19)  can  be  used  to  solve  for  the  new  gain 
relationship  with  modified  values  of  P  and  Jq,  represented  by  P*  and  Jq*. 


g=p*(J-Jo*). 


(4.20) 


This  is  done  by  solving  for  the  gain  as  a  function  of  distance  from  the  stripe 
center  and  current  density  through  the  stripe  contact,  Jtot. 

gain(x,Jtot)  =  P  (  (Jr(x)  Jtot  )  -  Jo)  •  (4.21) 


The  effective  gain  will  then  be  the  average  gain  over  the  region  where  there 
is  gain; 


gaineff(Jtot)  = 


1 

Xeff 


xeff 

J  gain(x,Jtot) 
0 


(I  X 


(4.22) 


where  Xeff  is  the  distance  from  the  center  of  the  stripe  to  the  point  where  the 
gain  becomes  zero.  Equation  (4.22)  makes  the  approximation  that  there  is  a 
uniform  density  of  photons  in  the  region  where  there  is  gain.  The  effective 
gain  can  be  found  for  other  values  of  Jtot  and  then  plotted.  The  slope  of  the 
curve  gives  the  new  value  of  P*,  the  gain  coefficient,  and  the  x-intercept  the 
new  transparency  current  density,  Jq*. 

Figure  4.2  shows  the  new  gain  curve  plotted  along  with  the  gain  curve 
for  no  current  diffusion  as  a  function  of  current  density.  This  is  for  the  case 
of  S  =  5  pm,  Ln  =  5  pm,  lo  =  1.2  pm,  P  =  4.4  cm/A,  and  Jq  =  145  A/cm^.  The  P  and  Jq 
values  are  taken  from  values  determined  for  the  SQW-60-300  laser.2 
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FIGURE  4.2  Gain  versus  current  density  through  the  stripe  contact. 


From  Figure  4.2  the  new  effective  gain  parameters,  |3*  and  Jq*. 
representing  the  slope  and  x-intercept  of  the  curve,  can  be  found.  Since  the 
graph  only  approximates  a  straight  line  these  values  must  be  interpolated.  It 

is  apparent  that  current  spreading  through  the  stripe  for  a  narrow  stripe 
width  has  dramatic  effects  on  the  effective  gain  of  the  laser.  Table  4.1  lists 
values  for  the  effective  gain  parameters  calculated  for  varying  spreading  and 
diffusion  lengths. 

TABLE  4.1 

Effective  Gain  Parameters  for  Various  Diffusion  Lengths  and  Stripe  Widths 


Diffusion 

Spreading 

Stripe 

P 

■lo 

Jo* 

Length 

Length 

Width 

(cm/Amp) 

(cm/Amp) 

(pm) 

(pm  ) 

(pm  ) 

5 

1.2 

5 

4.4 

145 

0.51 

225 

7 

1.2 

5 

*• 

'* 

0.415 

250 

7 

3 

5 

ti 

0.39 

325 

8 

5 

5 

0.32 

375 
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4.4  Line  shape  factor  and  the  gain  spectral  linewidth 
The  line  shape  factor  is  a  Lorentzian  equation  given  as 


Dv  = 


_ 1 _ 

l+((A.^-Xp{t))MXD]2 


(4.23) 


where 

AX.D  =  effective  gain  spectral  linewidth  parameter,  pm 
^p(t)  =  wavelength  at  the  peak  of  the  gain  spectrum  for  time  t 
=  wavelength  of  the  uth  laser  mode. 

The  spacing  between  adjacent  modes  in  a  cavity,  AX.^,  is  determined  by 


AXq  —  Xpo^/  (2ngL). 


(4.24) 


The  linewidth  parameter,  AXd.  has  the  effect  of  determining  the 
steady-state  relative  power  distribution  between  the  modes.  The  line  shape 
factor  equation  is  determined  by  fitting  the  Lorenztian  equation  to  the 
calculated  and  measured  spontaneous  emission  spectra  of  QW  GaAs  lasers.47-50 

4.5  Characterization  of  the  spontaneous  emission  factor 


Of  the  total  spontaneous  emission  emitted  in  all  directions,  only  a  small 
portion  becomes  amplified  by  the  laser  modes.  This  fraction  is  given  by  the 
spontaneous  emission  factor,  y,5  1,3 1,52, 5 3 


y= 


_ i^rXpo4 _ 

47i2nr^ngDdLAXD 


(4.25) 
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where  K  is  the  astigmatism  factor  and  Oj.  is  the  effective  refractive  index  of 
the  laser  medium.  The  factor  y  is  the  ratio  of  the  volume  occupied  by  one 
guided  mode  relative  to  the  whole  volume  receiving  photons  by  spontaneous 
emission. 

The  astigmatism  factor  is  used  to  account  for  changes  in  y  due  to  laser 
geometry.  As  will  be  discussed  in  Chapter  5,  K  is  affected  by  laser  length  and 
width.  It  also  changes  depending  on  whether  the  laser  is  gain-  or 
index-guided. 

The  spontaneous  emission  factor  affects  the  mode  power  distribution.  A 
larger  y  causes  a  larger  spread  of  the  power  distribution  of  the  modes.  The 
spontaneous  emission  factor  also  affects  the  time  it  takes  for  the  modal  power 
sharing  to  settle,  with  a  smaller  y  causing  the  modes  to  take  longer  to  stabilize 
in  power. 

4.6  Incorporation  of  temperature  effects 

The  peak  wavelength  of  the  gain  curve  is  dependent  on  the  current 
density^^  and  on  temperature  as  shown  in  the  following  equation, 

^i)-Ji.p(t)  =  u^l  -  5kpt(t)  +  AX  +  (u  •  AXc).  (4.26) 

In  this  equation,  AX  is  the  initial  detuning  of  the  gain  peak  from  the  dominant 
mode  (from  0-0.5  AXc),  and  u  is  a  constant  relating  the  carrier  density  to  the 
shift  in  the  gain  peak  (typically  106.4  nm)^^.  X^  and  Xp(t)  were  defined  earlier 

as  the  wavelength  of  the  uth  laser  mode  and  the  wavelength  at  the  peak  of  the 
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mode  distribution  for  time  t.  The  thermal  shift  of  the  peak  wavelength  of  the 
gain  curve,  6X.pt(t),  is  defined  as 

SXpt(t)  =  h  AT(t)  .  (4.27) 

where  h  is  the  experimental  shift  of  wavelength  per  ^C.  AT(t)  is  the 
temperature  shift  in  the  active  area  previously  described  by  equation  (3.1). 

The  gain  peak  shifts  due  to  carrier  density  changes  and  thermal  changes  are 
both  positive,  i.e.,  the  gain  peak  shifts  to  longer  wavelengths. 

The  carrier  density  term  becomes  zero  at  threshold  because  at  threshold 
and  at  higher  modulation  levels  the  carrier  density  stays  constant  except  for 
some  minor  oscillations.  Therefore,  the  factor  u  in  (4.26)  which  relates  the 
shift  in  the  gain  peak  due  to  carrier  density  changes  has  minimal  effect  in  the 
modeling  done  in  this  report  since  the  lasers  were  all  pulsed  from  threshold. 
The  carrier-induced  index  change  and  its  corresponding  wavelength  shift  is 
also  described  by  Manning,  Olshansky  and  Su.^'^ 

The  change  in  temperature  over  time,  AT(t),  is  determined  either 
experimentally  or  with  a  finite  element  analysis  of  the  thermal  transients  in 
the  active  area  of  the  laser  diode  during  a  pulse. To  determine  AT(t) 
experimentally,  the  modes  of  the  pulsed  laser  diode  are  separated  with  a 
spectrometer  and  the  shift  in  the  gain  curve  peak  during  the  pulse 
determined.  This  needs  to  be  done  for  each  separate  laser  at  each  power  level. 
A  plot  of  the  shift  in  the  gain  curve  with  time  can  be  made  and  the  data  fitted 
to  two  thermal  time  constants,  T1  and  T2,  where  T1  is  a  short  thermal  time 
constant,  and  T2  is  a  long  time  constant.  The  temperature  increase  in  the 
active  layer  for  time  constant  T1  is  kl  and  the  temperature  increase  in  the 
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active  layer  for  lime  constant  T2  is  k2.  Alternately  the  change  in  temperature 
over  time,  AT(t),  can  be  determined  with  a  finite  element  thermal  analysis  of 
the  diode  as  discussed  previously  in  Chapter  3. 

The  output  curves  from  the  rate  equation  analysis  in  this  paper  are 
compared  by  the  use  of  their  time  constants.  The  time  constant  of  the  curve 
comes  from  the  equation 

AT(t)  =  Tfinal  [l-e(-t^)]  .  (4.28) 

When  t  equals  x  the  value  is  at  63.2%  of  Tfinal.  its  final  value.  When  t  equals  2x 
the  value  is  at  86.5%  of  its  final  value. 

The  threshold  current  dependence  on  temperature  is  defined  as^ 


Iih(t)  =  Itho  elAT(t)/To]  .  (4.29) 

where  Itho  ‘s  the  threshold  current  at  the  initial  temperature,  lth(0  is  the 
threshold  current  at  a  time  t  later  as  the  temperature  changes  according  to 
Equation  (3.1),  and  Tq  is  the  active  region  characteristic  temperature  (-150K 
for  multi-quantum  well  diodes  and  -130K  for  single  quantum  well  diodes^-'*). 
The  threshold  current  increases  with  increasing  temperature.  For  example, 
for  To=150K  and  AT  =  5°K,  l^i  =  1.034Iiho.  For  AT  =  lO^K,  la,  =  1.0691iho- 

The  temperature  dependence  of  the  transparency  electron  current 
density,  Jq,  and  the  differential  gain  coefficient,  p,  can  be  derived  from  the 
threshold  current  dependence  using  (4.15).^^  If  the  temperature  dependence 
of  Jo  is  described  by  Jo(AT(t))=Joef^^(^)^oi  then  P(AT(t))= 

In  this  rate  equation  code,  all  the  spectral  shift  is  assumed  to  occur  via 
mode  hopping.  This  is  because  the  individual  modes  are  being  modeled.  The 
mode  chirping  can  easily  be  predicted,  though,  knowing  the  average  shift  in 
the  wavelength  of  the  individual  modes  with  temperature.  It  is  important  to 
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remember,  though,  when  comparing  the  modeling  to  the  experimental  data, 
that  if  mode  chirping  is  severe  enough,  the  experimental  data  will  not  give  an 
accurate  measure  of  the  mode  hopping,  because  the  mode  will  be  moving  out  of 
the  detector  area  during  the  pulse. 


4.7  Refractive  Indexes 

The  group  refractive  index,  ng,  is  defined  as 

,  ISr 

"g  -  "r  -  ^  dX  ’ 


(4.30) 


where  nr  is  the  effective  refractive  index  of  the  material.  The  group 
refractive  index  was  estimated  for  a  GaAlAs  SQW  laser  at  an  energy  of  1.455  eV 
to  be  3.87.^^  Larsson,  et.al.,56  report  values  for  the  group  index  as  3.45  to  3.8 
for  GRINSCH  SQW  lasers  with  cavity  lengths  ranging  from  300  to  600  pm.  This 
value  concurs  with  values  obtained  from  the  measurements  of  the  separation 
of  the  laser  diode  modes  with  a  spectrometer  and  equation  (4.24).  These  values 
are  listed  in  Table  5.2  of  Chapter  5  for  several  of  the  lasers  tested  in  this  report. 
The  effective  refractive  index  of  the  GaAs  laser  medium  for  laser  energy  at 
1.5  eV  is  found  to  be  near  3.65. ^  2 


4.8  Confinement  factor 

The  confinement  factor,  f,  is  the  fraction  of  the  mode  energy  contained 
in  the  active  region.  It  also  is  an  approximate  ratio  of  the  electron  volume  in 
the  active  area  to  the  volume  occupied  by  the  optical  modes,  and  is  defined  as 
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r  = 


dll 


h 


E 
-dll 


(4.31) 


where  Ey  is  the  electric  field  and  d  is  the  active  area  thickness.^  The 
confinement  factor  can  be  approximated  by  the  of  the  effective  optical 
width,  s, 


s  = 


\eI  d  X 


'max 


(4.32) 


where  E^iax  *he  electric  field  at  the  center  of  the  active  area.  The  effective 
optical  width  is  close  to  the  full  width/half  maximum  of  the  optical  power 
profile.  The  confinement  factor  can  then  be  approximated  by 


r.i 

S 


(4.33) 


A  simple  relationship  calculating  f  from  the  refractive  indexes  of  the 
material,  the  active  layer  thickness,  and  the  wavelength  can  be  found  in 
Agrawal  and  Dutta.34  This  relationship  is  an  approximation  assuming  most  of 
the  power  is  in  the  fundamental  transverse  mode.  The  normalized  waveguide 
thickness  is  now  introduced; 


D  =  ko  (Mg  -  uj)-^  d 


(4.34) 
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where  and  are  the  refractive  indices  of  the  active  layer  and  the 

a  C 

cladding  layers,  k©  is  the  wave  number  of  the  mode. 


2  71 
^0 


(4.35) 


and  X-o  is  the  free  space  wavelength.  The  confinement  factor  can  then  be 
approximated  as34 


r  =  d2/(2  +  d2) 


(4.36) 


for  D  in  pm.  For  quantum  well  lasers,  where  D«  1,  this  equation  can  be  further 
simplified  to  F  =  D^/2  or 


r  =  2  tc2  (p2  .  d^/xj.  (4.37) 

Table  I  shows  some  calculations  of  the  confinement  factor  and  effective 
optical  thicknesses  for  various  refractive  indexes,  wavelengths,  and  active 
layer  thicknesses. 


Table  4.2 

Calculated  Confinement  Factors  and  Effective  Optical  Thicknesses 


Xq.  pm 

d,  pm 

r 

s,  pm 

0.8 

1.0 

3.51 

3.22 

0.98 

1.02 

0.8 

0.1 

3.51 

3.22 

0.376 

.266 

0.8 

.06 

3.51 

3.22 

.178 

.337 

0.8 

.02 

3.51 

3.22 

.0235 

.851 

0.8 

.01 

3.51 

3.22 

.006 

1.67 
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Kressel  and  Butler^  ^  and  Thompson^  get  similar  results  as  those  shown  in 
Table  I.  As  can  be  seen,  the  effective  optical  thickness,  s,  increases 
dramatically  for  active  layer  sizes  less  than  0.1  pm. 

Quantum  well  structures  allow  high  gain  at  low  carrier  densities  but 
have  the  disadvantage  of  a  small  optical  confinement  factor.^^  Two  means  of 
combating  this  are  the  use  of  multi-quantum  wells  and  a  graded  index 
confinement  structure. 

The  optical  confinement  factor  for  graded-index  structures  must  be 
calculated  numerically.  Wagner  et.al.^  calculates  this  for  various  laser 
lengths  and  widths.  For  a  laser  with  a  quantum  well  active  layer  thickness  of 
10  nm,  the  confinement  factor  was  0.032,  giving  an  effective  optical 
thickness,  s,  of  0.312  pm.  This  is  a  considerable  improvement  over  the 
non-graded  structure.  For  the  0.014  pm  active  layer  thickness  SQW  GRINSCH 
lasers  used  in  this  report,  a  confinement  factor  of  0.045  was  used  which  gives 
the  same  effective  optical  thickness  of  0.312  pm. 

Agrawal  and  Dutta^^  calculate  the  confinement  factor  for 
multi-quantum  -well  lasers  as 


r  =  Y 


Nada 

N  ada  +  N^db 


(4.38) 


where 


Y  =  2  rt^fNada  +  Nbdb)^ 


Nadalta  NbdbPb 
Nada  +  Nbdb 


and  Na  and  Nb  are  the  number  of  active  and  barrier  layers.  Using  these 

formulas  F  is  calculated  to  be  0.147  for  a  0.808  pm  wavelength  laser  with  four 
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active  layers  each  13  nm  thick  and  three  barrier  layers  each  4  nm  thick.  In 
this  case  s  will  be  0.436  trm . 


4.9  Cavity  loss  coefficient 

The  effect  of  the  loss  coefficient,  Oi,  in  the  late  equations  is  to  change 
the  overall  output  power  of  all  the  modes.  Raising  ai  will  lower  the  output 
power.  The  main  cavity  loss  mechanisms  are  free-carrier  absorption  in  the 

quantum  well  and  in  the  graded  and  confining  layers.  The  intrinsic  loss 
coefficient  is  related  to  these  by  the  confinement  factor, 


tti  =  ttqw  r  +  ttcl  (1  -  f). 


(4.39) 


The  free  carrier  absorption  coefficients  in  the  quantum  well  and  graded 
and  confining  layers,  aqw  and  Ud.  have  values  in  the  range  of  30  cm'^  and 
2  cm*^.  In  Wagner,  Waters,  et.al.,^  values  for  oj  are  found  for  a  SQW, 
GRINSCH  structure  laser.  They  obtain  values  from  2.8  to  5.5  cm*^  for  different 
structures. 

The  external  quantum  efficiency,  qe-  can  be  related  to  the  internal 
quantum  efficiency,  qi,  by  the  equation: 


ln(l/R) 

+  ln{l/R) 


(4.40) 


The  internal  quantum  efficiency  of  these  lasers  is  known  to  be  near 
1.0. '2  xhe  differential  quantum  efficiency  is  the  slope  of  the  optical  power 
versus  injection  current  graph  (in  W/A)  divided  by  the  bandgap  energy  in  eV 
where  the  bandgap  energy  in  eV  is  equal  to  1.24Apo.  with  X.po  in  microns. 
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4.10 


Spontaneous  Electron  Lifetime 


The  carrier  (or  electron)  lifetime,  Te.  is  the  carrier  lifetime  in  the  active 
region.  It  is  the  inverse  of  the  carrier  recombination  rate  and  is  defined  by 

—  =  Anr  +  Biie  +  Cne^ 

where 

Anr  =  non-radiative  recombination  rate  due  to  traps  or  surface 
recombination,  s‘^* 

B  =  radiative  recombination  coefficient,  cm^/s, 

C  =  Auger  recombination  coefficient,  cm^/s  . 


The  carrier  lifetime  is  related  to  the  carrier  lifetime  due  to  radiative 
recombination  processes,  tr,  and  the  carrier  lifetime  due  to  nonradiative 
recombination  processes,  Xnr.  by. 


te 


(4.42) 


where 


1  1  -7 

—  =  B  ne  and  =  Anr  +  C  ne"^. 
tr  tnr 


Olshansky,  et.  al.,  have  reported  values  for  the  Anr.  B.  and  C  coefficients 
for  AlGaAs  LEDs  and  InGaAsP  lasers.^^  They  found  that  the  Auger  coefficient 
is  negligible  for  AlGaAs  LEDs.  It  should  be  noted  that  the  carrier  lifetime  will 
remain  a  constant  once  threshold  is  reached  because  at  that  point  and  at 
higher  current  levels  the  electron  density,  ne,  remains  a  constant.  The  rate 
equation  modeling  in  this  report  is  d,;signed  to  model  lasers  that  are  pulsed 
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from  threshold  and  so  a  single  value  for  the  carrier  lifetime  at  threshold  can 
be  used.  If  it  was  desired  to  model  the  lasers  pulsed  from  zero  current  the  Anr 
and  B  coefficients  would  need  to  be  used.  The  Anr  and  B  coefficients  can  be 
found  by  measurements  of  the  carrier  lifetime  versus  current  as  described  by 
Olshansky.  Since  the  Anr  coefficient  is  usually  on  the  order  of  5  to  10  times 
smaller  than  the  Bne  term  at  threshold,  tr  is  approximately  equal  to  Xe  and 
there  is  only  a  relatively  small  discrepancy  if  Xe  is  used  for  Xr  in  Equation  (4.2) 
so  that  the  equations  can  be  normalized  more  easily. 

The  carrier  lifetime  for  the  MQW -4-250  laser  can  be  found 
experimentally  from  the  fall  time  of  a  pulse.  The  diode  was  pulsed  below 
threshold  and  the  fall  time  measured.  For  this  work  a  Berkeley  Nucleonics 
Model  6040  Pulse  Generator  with  a  Model  20 IE  Module  with  150  ps  rise  and  fall 
times  was  used  to  pulse  the  laser.  An  Antel  Optronics  Model  AR-S5  APD  with  a 
<90  ps  rise  time  was  used  to  detect  the  signal.  For  the  MQW-4-250  single-stripe 
laser  the  fall  time  was  found  to  be  approximately  1.2  ns. 

The  carrier  lifetime  can  also  be  calculated  from  the  expression^  ^ 

Xd  =  to  +  Xe  ln[I/(I  -  Iih)]  .  (4.43) 

where  xd  is  the  delay  time  between  the  onset  of  the  current  pulse  and  the 
output  from  the  laser,  and  to  is  a  fixed  arbitrary  delay  caused  by  cables.  The 
pulse  current  is  measured  from  an  initially  unbiased  value.  With  this 
equation,  if  several  values  are  taken  and  plotted,  the  electron  lifetime  can  be 
found  from  the  slope  of  the  curve  of  (xd)  versus  (xe  ln[I/(I  -  Ith)])-  Elliott, 
et.al.,''  calculated  the  lifetime  this  way  for  the  10-siripe  MQW-100-250  laser 
and  found  a  value  of  1.9  ns.  Streifer.  et.al.,  use.s  ilie  value  of  1.6  ns  for  a  MOW 
laser  at  threshold.^  ® 
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Changing  the  carrier  lifetime  in  the  rate  equations  program  changes 
the  oscillation  frequency.  This  is  because  it  is  directly  related  to  the  stimulated 
emission  factor  as  shown  in  Equation  (4.12).  A  larger  value  of  te  decreases  the 
number  of  oscillations  per  nanosecond. 

4.11  Limitations  of  the  rate  equation  model 

Several  factors  limit  the  accuracy  of  our  model. 

Firstly,  the  gain/current  relationship  is  assumed  to  be  linear  when  in 
fact  there  are  non-linearities  present. 

Secondly,  the  model  does  not  take  into  account  mode  chirping.  It 

assumes  the  modes  stay  at  a  constant  wavelength  and  only  the  power 
distribution  between  modes  changes. 

Values  for  y,  Xsp*  and  a  are  approximate  and  will  vary  even  among 
lasers  of  the  same  type. 

However,  with  all  these  inaccuracies,  the  rate  equation  model  is  still 
useful  as  a  tool  to  obtain  rough  measurements  of  the  behavior  of  pulsed  laser 
diodes,  particularly  when  various  design  parameters  need  to  be  compared. 

A  complete  listing  of  the  rate  equation  computer  program  in 
Mathematica  for  a  Macintosh  computer  is  given  in  Appendix  A. 
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CHAPTER  5 


USING  THE  RATE  EQUATION  MODEL 

5.1  Introduction 

The  rate  equation  model  can  be  used  with  or  without  including  thermal 
effects.  Without  thermal  effects  the  model  can  be  run  to  get  information  on 
relaxation  oscillation  frequency,  modal  power  distribution,  total  output  power, 
and  the  settling  time  of  the  modal  power  fluctuations  due  to  mode  buildup-time 
transients. 

When  temperature  dependence  is  added,  the  relaxation  oscillation 
frequency  and  total  output  power  remain  the  same.  The  change  is  that  now 
the  peak  of  the  gain  curve  will  shift  and  the  modes  will  shift  in  power 
accordingly.  The  longer  wavelength  modes  will  gain  power  and  the  shorter 
wavelength  modes  will  lose  power. 

5.2  Modeling  the  laser  diodes  without  incorporating  thermal 
effects 

Table  5.1  lists  the  input  parameters  for  the  lasers  used  in  the 
simulations.  Wavelength  and  spacing  between  the  individual  modes,  AXc, 
taken  with  an  LS-2  spectrometer,  are  also  presented.  These  were  used  to 
calculate  values  for  the  group  index,  ng.  A  value  of  1.8  ns''  was  used  for  Xgp, 
0.025  |im  for  AX,d.  and  0.3  and  0.95  for  R1  and  R2.  Section  4.8  describes  how  the 
confinement  factors,  f,  were  obtained.  Section  4.3  describes  how  the 
differential  gain  coefficient,  p,  and  the  transparency  current  density,  Jo,  were 
obtained. 
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TABLE  5.1 


Individual  Input  Parameters  for  the  Laser  Diodes 


Parameter 

SOW-60-300 

SOW-60-600 

D  (urn) 

60 

60 

5 

4 

L  (urn) 

300 

600 

300 

250 

d  (pm) 

0.014 

0.014 

0.014 

0.064 

r 

0.0452 

0.0452 

0.0452 

0.147 

B  (cm/A) 

4.4 

4.4 

0.51 

0.12 

Jo  (A/cm^) 

145 

145 

225 

565 

"r 

3.65*2 

3.65*2 

3.65*2 

3.6 

^00  (uml 

0.814 

0.822 

0.819 

0.80 

Ith  (mA) 

100 

100 

75 

49 

Differential  Quantum 
Efficiency 

0.5 

0.623 

0.42 

^Do  (um)  (from  LS-2) 

0.812 

0.818 

0.811 

0.802 

AXc  fnml  (from  LS-2'i 

0.312 

0.137 

0.305 

0.30 

iLfi _ 

3.52 

3.96 

3.59 

4.3 

Simulations  of  the  modal  transient  response  without  thermal  effects  for 
four  lasers  were  performed.  Table  5.2  shows  the  results  along  with  values 
obtained  experimentally.  The  relative  output  power  of  the  individual  modes 
for  the  lasers  was  not  determined  experimentally  but  can  be  compared 
qualitatively  from  the  experimental  data  presented  in  Figures  2.4  through  2.11. 
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TABLE  5.2 


Rate  Equation  Program  Results  for  the  Laser  Diodes 


Parameter 

Calc. 

Value 

Exper 

Value 

Calc. 

Value 

Calc. 

Value 

Exper 

Value 

Calc. 

Value 

Laser 

SOW-60-300 

SOW-60-600 

SOW-5-300 

BX7!CKRE?!I 

Current,  rr.A 

200 

200 

1  28 

80 

Intrinsic  Loss,  oj,  cm"  ^ 

29 

14 

39 

60 

Astigmatism  Factor,  K 

68 

32 

2.5 

2 

Total  Output  Power,  mW 

80 

83 

65 

68 

46 

46 

20 

19 

Threshold  Current,  mA 

100 

74 

1  00 

99 

75 

44 

49 

38 

Oscill.  Frequency,  ns‘  ^ 

3.2 

4.0 

2.4 

2.4 

1 .7 

3.8 

warn 

3.3 

Mode  buildup  time, 
center  mode,  ns 

-10 

10 

-30 

30 

-10 

1  5 

-8 

15 

Output  Power,  mW: 

Mode  1 

- 

23 

1  2 

- 

23 

- 

9 

Mode  2 

- 

1  3 

. 

1  0 

7 

3.8 

Mode  3 

- 

7 

- 

6 

1  .5 

. 

1 

Mode  4 

- 

3 

- 

4 

- 

0.5 

- 

0.5 

When  the  rate  equation  runs  are  being  done,  the  gain  parameters, 
carrier  lifetime,  and  laser  dimensions  are  entered  into  the  program  and  then 
trial  and  error  runs  are  made  to  fit  the  calculated  results  to  the  experimental 
data  by  adjusting  the  astigmatism  factor  and  the  intrinsic  loss  coefficient.  The 

astigmatism  factor  (as  part  of  the  spontaneous  emissitn  factor,  y)  affects  the 
mode  power  distribution  and  to  a  small  extent  the  oscillation  frequency  of  the 
relaxation  oscillations.  A  larger  7  causes  a  larger  spread  of  the  power 
distribution  of  the  modes  and  a  slightly  higher  oscillation  frequency.  The 
spontaneous  emission  factor  also  affects  the  lime  it  takes  for  the  modal  power 
sharing  to  settle,  with  a  smaller  7  causing  the  modes  to  lake  longer  to  stabilize 
in  power.  The  effect  of  the  loss  coefficient,  ai,  in  the  rate  equations  is  to 
change  the  overall  output  power  of  all  the  modes.  Raising  Oi  will  lower  the 
total  output  power.  These  values  of  the  astigmatism  factor  and  the  intrinsic 
loss  coefficient  are  only  approximations,  though,  due  to  the  approximate 
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nature  of  several  of  the  other  input  parameters,  such  as  the  confinement 
factor,  the  electron  lifetime,  and  the  gain  coefficients. 

Close  agreement  is  obtained  between  the  experimental  and  calculated 
values  for  the  SQW-60-300  laser  for  the  total  output  power,  the  relaxation 
oscillation  frequency,  and  the  mode  buildup  time.  Figure  5.1a  shows  the  output 
from  four  of  the  dominant  modes  for  the  SQW-60-300  laser.  This  can  be 
compared  to  the  experimental  data  in  Figure  2.4.  Mode  0  is  the  center  mode, 
and  modes  1  through  3  are  side  modes  which  have  identical  counterparts  on 
the  other  wavelength  side  of  mode  0.  The  mode  buildup-time  represents  the 
time  it  takes  for  the  center  mode  to  reach  86%  (two  time  constants)  of  the  final 
value.  The  other  modes  settle  faster,  at  rates  that  are  proportional  to  their 
distance  from  the  center  mode.  Figure  5.1b  shows  the  electron  concentration 
as  a  function  of  time  during  the  pulse.  All  the  rate  equation  runs  are  started 
from  a  current  slightly  below  threshold.  As  can  be  seen,  the  electron 
concentration  rises  to  the  threshold  value  and  then  oscillates  around  this 
value  before  leveling  out. 

Figure  5.2  shows  the  output  from  six  of  the  dominant  modes  for  the 
SQW-60-600  laser  diode.  A  smaller  astigmatism  factor,  K,  for  the  spontaneous 
emission  factor,  y,  was  needed  for  this  laser  to  match  the  experimental  data. 

The  smaller  y  needed  for  this  laser  is  related  to  the  fact  that  with  the  longer 
length,  the  allowed  angular  spread  of  the  individual  modes  is  less  so  that  a 
smaller  percentage  of  the  spontaneous  emission  (which  is  emitted  in  all 
directions)  will  go  into  the  lasing  modes. 

The  mode  buildup-time  for  the  SQW-60-600  laser  is  considerably  longer 
than  for  the  SQW-60-300  laser.  This  is  related  to  the  fact  that  it  has  a  smaller 
spontaneous  emission  factor  and  also  to  the  fact  that  the  photon  lifetime  for 
this  laser  will  be  longer  because  of  the  longer  length  and  so  the  competition 
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FIGURE  5.1  Rate  equation  modeling  for  the  SQW-60-300  laser  diode 
showing  (a)  output  power;  and  (b)  relative  electron  concentration;  from 
four  lasing  modes. 
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FIGURE  5.2  Rate  equation  modeling  for  the  SQW-60-600  laser  diode 
showing  output  from  six  lasing  modes. 
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between  the  modes  takes  longer  to  resolve.  The  mode  spacing  for  this  laser  is 
half  that  of  the  SQW-60-300  laser  and  the  output  power  is  spread  out  more 
between  the  modes. 

Figure  5.3  shows  the  output  from  the  four  dominant  modes  for  the 
SQW-5-300  laser  diode  using  the  effective  values  of  P  and  Jo  of  0.51  cm/A  and 
225  A/cm^  calculated  using  a  diffusion  length  of  5  pm  and  current  spreading 
length  of  1.2  pm.  The  total  output  power  and  mode  buildup  time  matched  the 
experimental  data  well,  however,  the  calculated  oscillation  frequency  for  this 
laser  of  3.8  GHz  did  not  match  the  experimental  value  of  1.7  GHz  well.  This 
indicates  that  either  the  effective  value  of  p  is  too  high  or  that  the  carrier 
lifetime,  Xe,  is  too  low.  Using  effective  values  of  P  and  Jo  of  0.32  and  375, 
determined  from  the  current  spreading  calculation  with  a  diffusion  length  of 
8  pm  and  current  spreading  length  of  5  pm,  a  value  of  2.6  GHz  was  obtained  for 
the  oscillation  frequency.  With  these  longer  diffusion  and  current  spreading 
lengths,  a  calculated  current  threshold  of  66  mA  was  obtained  which  is  closer 
to  the  experimental  value  of  75  mA  than  the  value  of  44  mA  obtained  using  the 
calculated  diffusion  and  current  spreading  lengths.  Using  these  longer 
diffusion  and  current  spreading  lengths  is  not  reasonable,  however,  since  the 
experimental  data  taken  by  imaging  the  laser  facet  below  threshold  does  not 
indicate  that  there  is  that  much  spread  in  the  current  profile.  It  is  apparent 
that  some  other  factor,  such  as  a  change  in  the  carrier  lifetime,  is  playing  a 
role  in  the  behavior  of  these  narrow  stripe  lasers. 

The  SQW-5-300  laser  needed  a  much  smaller  astigmatism  factor,  K,  of  2.5 
to  match  the  steady  state  distribution  of  power  between  the  modes.  This  is  due 
to  the  fact  that  since  it  is  a  narrow  stripe,  the  allowed  angular  spread  of  the 
individual  modes  is  less  than  for  a  broad  stripe  laser  and  a  smaller  percentage 
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FIGURE  5.3  Rate  equation  modeling  for  the  SQW-5*300  laser  diode 
showing  output  from  four  lasing  modes. 
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of  the  spontaneous  emission  (which  is  emitted  in  all  directions)  will  go  into 
the  allowed  lasing  modes. 

Table  5.1  also  lists  the  input  parameters  used  to  model  the  MQW-4-250 
narrow  stripe  laser  diode.  A  MQW  laser  typically  consists  of  four  0.013  pm 
thick  layers  of  Ga.94 A  1.06^5  separated  by  three  0.004  pm  thick  layers  of 
Ga.8Al.2As. ^ ®  An  effective  thickness  of  0.064  pm  and  confinement  factor  of 
0.147  were  used  to  model  this  laser. 

The  experimental  relaxation  oscillations  from  the  MQW-4-250  laser  were 
det  ;rmined  in  a  paper  by  Elliott  and  DeFreez.^^  The  relaxation  oscillations  had 
a  frequency  of  1.7  GHz  (1.7  oscillations/ns)  for  a  current  1.3  times  the 
threshold  current  and  4.0  GHz  (4  oscillations/ns)  for  a  current  twice 
threshold.  This  data  can  be  extrapolated  for  a  current  1.6  times  threshold  to 
give  2.7  oscillations/ns. 

The  rate  equation  results  for  the  MQW-4-250  are  presented  in 
Table  5.2.  Effective  values  of  P  and  Jq  of  0.2  and  750  were  found  to  give  results 
cli/se  to  those  found  experimentally.  (Intrinsic  values  of  P  and  Jq  could  not  be 
fc  md  for  this  laser.)  Figure  5.4  shows  the  output  from  the  four  dominant 
m  ides.  The  value  of  2.0  for  K  needed  for  this  laser  is  close  to  the  value  of  3 
calculated  in  a  report  by  Hausser,  et.  al.,^^  for  narrow  stripe  MQW  lasers. 

All  the  runs  were  done  using  19  modes  in  the  calculation.  When  only  9 
modes  were  used  the  results  were  similar.  With  9  modes  the  oscillations  were 
more  prominent  but  the  modes  settled  to  approximately  the  same  powers  and 
the  average  power  was  approximately  the  same. 

A  run  was  also  done  for  the  SQW-60-300  laser  with  the  gain  curve 
peaked  between  two  modes.  In  this  case  the  results  are  very  similar  except  for 
the  fact  that  there  are  now  two  identical  center  modes  with  slightly  less  power 
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FIGURE  5.4  Rate  equation  modeling  for  the  MQW -4-250  laser  diode 
with  20  mW  output  power  and  a  20  ns  pulse. 
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than  the  one  center  mode  for  the  case  of  the  gain  peak  centered  on  one  mode. 
The  settling  times  (or  mode  buildup  time)  for  the  modes  remained 
approximately  the  same. 

Although  the  gain  saturation  factor,  e.  was  not  used  in  the  rate  equation 
runs,  values  for  e  could  be  obtained  by  trial  and  error  that  damp  out  the 
oscillations  in  about  two  nanoseconds.  A  value  of  2.0x10'^  pm  ^  was  obtained 
for  the  SQW-5-300  laser.  24x10-6  pm  3  for  the  SQW-60-300  and  the  SQW-60-600 
lasers,  and  3.5x10-6  pm^  for  the  MQW-4-250  laser. 

5.3  Modeling  the  laser  diodes  with  thermal  effects 

As  stated  earlier,  when  temperature  dependence  is  added  to  the  rate 
equation  modeling,  the  relaxation  oscillation  frequency  and  total  output  power 
remain  the  same.  The  change  is  that  now  the  peak  of  the  gain  curve  will  shift 
with  time  and  the  modes  will  shift  in  power  accordingly.  The  longer 
wavelength  modes  will  gain  power,  and  the  shorter  wavelength  modes  will 
lose  power. 

Table  3-8  lists  the  thermal  time  constants  for  the  SQW-60-300  and 
SQW-5-300  lasers  derived  from  the  finite  element  analysis.  The  parameters  kl, 
k2,  Tl,  and  T2  are  dependent  on  how  much  drive  current  is  applied  to  the  diode 
and  so  must  be  determined  for  a  specific  value  of  current. 

The  characteristic  temperature  used  was  125  °K  2.  A  value  of 
0.25  nm/°K  was  used  for  h. 

Figures  5.5  and  5.6  show  the  rate  equation  modeling  with  thermal 
effects  for  the  SQW-60-300  laser  for  80  mW  power  output  (100  mA  modulation 
current)  and  80  ns  and  1000  ns  pulses  using  the  finite  element  analysis 
derived  thermal  parameters.  These  can  be  compared  to  the  experimental  data 
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FIGURE  5.5  Rate  equation  modeling  for  the  SQW-60-300  laser  with  thermal 
effects  derived  from  the  finite  element  analysis  for  80  mW  output  power  and 
an  80  ns  pulse. 
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in  Figures  2.4  and  2.5.  This  laser  had  fairly  small  thermal  fluxes  for  the  80  mW 
power  level  and  this  is  confirmed  by  the  modeling.  The  thermal  transients  for 
the  80  ns  pulse  are  too  small  to  be  detected.  Small  thermal  transients  of  the 
same  magnitude  are  observable,  though,  for  the  1000  ns  pulse  in  the  modeling 
and  the  experimental  data. 

Figure  5.7  shows  the  rate  equation  modeling  for  the  SQW-5-300  laser  for 
46  mW  power  output  and  an  80  ns  pulse  using  the  finite  element  analysis 
derived  thermal  parameters  with  a  FWHM  of  5  pm  for  the  current  profile.  As 
can  be  seen  by  comparing  this  data  to  the  experimental  data  in  Figure  2.8,  the 
experimental  thermal  transients  are  reasonably  close  to  the  calculated  ones. 
Figure  5.8  shows  the  SQW-5-300  laser  modeled  with  46  mW  optical  power  output 
and  a  1000  ns  pulse.  Once  again,  the  experimental  thermal  transients  look 

reasonably  close  to  the  calculated  ones,  though,  a  direct  comparison  cannot  be 
made  to  the  experimental  data  of  Figure  2.9  because,  as  stated  earlier,  the 
transients  for  this  laser  are  of  a  high  enough  magnitude  to  cause  the  modes  to 
shift  through  the  detector  area  during  the  pulse,  making  the  degree  of  mode 
hopping  appear  larger. 

5.4  General  discussion  of  rate  equation  modeling  results 

The  rate  equation  results  agree  with  the  experimental  data  in 
describing  the  two  types  of  transients  occurring  in  these  gain-guided  lasers. 

Mode-buildup  transients  are  more  severe  and  occur  at  the  beginning  of  the 
pulse.  Their  duration  is  dependent  on  the  laser  length,  with  15  ns  transients 
occurring  in  300  pm  length  lasers,  and  30  ns  transients  occurring  in  600  pm 
length  lasers. 
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FIGURE  5.7  Rate  equation  modeling  for  the  SQW-5-300  laser  with  thermal 
parameters  derived  from  the  finite  element  analysis  for  46  mW  output 
power,  an  80  ns  pulse,  and  a  FWHM  of  5  pm  for  the  current  profile. 
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FIGURE  5.8  Rale  equation  modeling  for  the  SQW-5-300  laser  with  thermal 
parameters  derived  from  the  finite  element  analysis  for  46  mW  output 
power,  a  1000  ns  pulse,  and  a  FWHM  of  5  pm  for  the  current  profile. 


The  rate  equation  program  was  found  to  give  a  good  match  to  the 
experimental  data  for  the  broad  stripe  lasers.  The  minor  discrepancies  can  be 
attributed  to  the  inherent  inaccuracy  of  some  of  the  simplifying  assumptions 
and  the  approximate  nature  of  several  of  the  input  parameters.  The  thermal 
transients  obtained  from  the  rate  equation  analysis  matched  the  experimental 
data  well.  The  magnitude  of  the  calculated  thermal  transients  are  dependent 
not  only  on  the  finite  element  analysis  derived  time  constants,  but  also  on  the 
spontaneous  emission  factor.  The  spontaneous  emission  factor  has  a  direct 
effect  on  the  relative  power  distribution  between  the  modes.  If  the  power  is 
spread  more  evenly  between  several  modes  (using  a  larger  spontaneous 
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because  a  shift  in  the  g.^in  peak  does  not  cause  as  great  a  shift  in  the 
magnitude  of  the  modes. 

For  the  narrow  stripe  lasers  the  match  of  the  experimental  data  to 
the  calculated  data  was  not  as  good.  Longer  diffusion  and  current 
spreading  lengths  than  those  calculated  were  needed  to  adjust  the  gain 
parameters  to  match  the  experimental  relaxation  oscillation  frequency  and 
threshold  current.  However,  these  longer  lengths  were  not  confirmed  by 
experimental  data  taken  by  imaging  the  laser  facet  below  threshold.  The 
experimental  data  indicated  that  the  current  profile  narrows  considerably 
near  threshold.  A  reasonable  match  of  the  experimental  and  calculated 
thermal  transients  was  obtained  by  using  a  FWHM  of  5  pm  for  the  current 
profile.  Modeling  of  narrow  stripe  lasers  is  an  area  for  further 


investigation. 
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CHAPTER  6 


CONCLUSIONS 

In  this  report  the  axial  mode  instabilities  for  several  laser  diode 
structures  were  examined.  In  Chapter  2,  experimental  data  was  taken  and  the 
general  results  showed  two  types  of  modal  transients  present.  Mode  buildup 
transients  occur  over  a  relatively  short  period  of  time  (10-60  ns)  that  is 
dependent  on  the  photon  lifetime  (and  consequently  laser  length,  with 
shorter  lengths  having  smaller  transients).  Thermal  transients  occur  over  a 
period  of  several  hundred  nanoseconds. 

The  broad-stripe  GRINSCH  SQW  laser  diodes  had  negligible  thermal 
transients  for  the  peak  modulation  depth  of  100  mA.  The  MQW  narrow  stripe 
laser  and  SQW  narrow-stripe  laser  had  large  thermal  transients. 

In  Chapter  3,  the  steady  state  temperature  increase  of  the  laser  diodes 
and  the  temperature  changes  in  the  active  area  during  a  current  pulse  were 
modeled  with  a  finite  element  code.  A  detailed  analysis  of  the  laser  was  done 
2-dimensionally,  with  a  3-dimensional  analysis  performed  to  determine  the 
thermal  impedance  of  the  heal  sink. 

In  Chapter  4,  the  thermal  modeling  was  combined  with  rate  equations  to 
give  a  complete  model  of  how  the  laser  modes  behave  during  modulation.  The 
rate  equation  model  without  thermal  effects  gave  information  on  relaxation 
oscillation  frequency,  threshold  current,  modal  power  distribution,  total 
output  power,  and  the  settling  time  of  the  modes  due  to  mode  buildup 
transients.  Adding  thermal  effects  shifted  the  power  distribution  of  *he  modes 
according  to  the  entered  thermal  lime  constants  in  a  very  predictable  fashion. 

The  rate  equation  results  agreed  with  the  experimental  data  in 
describing  the  two  types  of  transients  occurring  in  these  gain-guided  lasers. 
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Mode  buildup-time  transients  are  more  severe  and  occur  at  the  beginning  of 
the  pulse.  Their  duration  is  dependent  on  the  laser  lengths,  with  10-30  ns 
transients  occurring  in  300  pm  length  lasers,  and  20-60  ns  transients 
occurring  in  600  pm  length  lasers.  The  rate  equation  program  was  found  to 
give  a  good  match  to  the  experimental  data  for  the  broad  stripe  lasers.  The 
minor  discrepancies  can  be  attributed  to  the  inherent  inaccuracy  of  some  of 
the  simplifying  assumptions  and  the  approximate  nature  of  several  of  the 
input  parameters.  Using  the  rate  equation  model  with  the  finite  element  code 
derived  thermal  parameters  was  found  to  give  a  good  match  to  the 
experimental  data. 

A  method  was  found  to  solve  for  the  gain  constant  and  the  transparency 
current  density  for  narrow-stripe  lasers.  Current  spreading  in  the 
narrow-stripe  lasers  was  found  to  have  a  large  effect  on  these  constants. 
Diffusion  and  current  spreading  lengths  longer  than  those  calculated, 
however,  were  necessary  to  match  the  experimental  data  on  the  relaxation 
oscillation  frequency  and  the  threshold  current.  This  was  not  consistent  with 
the  thermal  modeling  and  the  experimental  data  taken  by  imaging  the  laser 
facet,  which  indicates  that  the  current  spreading  region  should  be  less  than 
that  calculated.  A  good  match  to  experimental  thermal  transients,  however, 
was  obtained  using  the  thermal  modeling  with  the  rate  equations. 

Further  works  needs  to  be  done  on  modeling  the  narrow  stripe  lasers, 
however,  it  should  be  noted  that  these  laterally  gain-guided  narrow  stripe 
lasers  are  not  used  as  frequently  nowadays,  because  index-guided  or  buried 
heterostructure  lasers,  which  have  lateral  current  confinement,  are  more 
commonly  available  and  arc  preferred  over  gain-guided  devices  particularly 
for  communication  purposes. 
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A  broad-stripe  GRINSCH  SQW  laser  diode  with  a  200-250  pm  length  would 
be  a  good  structure  for  minimizing  both  mode  buildup-time  transients  and 
thermally-induced  transients.  This  type  of  structure  could  output  100  mW  of 
output  power  with  minimal  modal  instability.  Going  to  higher  output  powers 
could  still  be  a  problem,  though,  since  the  thermal  transients  are  directly 
related  to  modulation  depth. 
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(• 

T2-1500 

(• 

E1-.29S 

(• 

E2-.41 

(• 

TO-125 

(• 

(• 

jon-0 

(• 

<• 

(*  fPSCTRAL  9ARJ 

Rnndtypa-1 

(• 

R0-.02S 

(• 

Rl-.02< 

(• 

R2-.014 

(• 

lnmbdap-0.912 

(• 

v-.COO 

(• 

(• 

(• 

1^9 

(• 

SU-. 01044 

(• 

u-0 


(*  valoeity  of  light,  aierona  par  na  *) 

(*  valoeity  of  light,  aierona  par  aao  *) 

(*  alactron  eharga  C/aloctron  *) 

0*9  (*  flanck*#  oonatant,  Joula  aa  *) 

7t$  tARMaTEM  •) 

gain  apactral  ahift  for  t«aip  in  aierena/E  .0003*) 
aat  Ton  to  1  to  ineluda  taa^ratura  affacta,  *) 
aat  to  0  to  ignora  affacta  *) 

taap.  inc.in  activa  layer  froa  Xtb  to  Xp  at,  na  *) 

Tbaraal  tiaa  oonatant  of  activa  layer  in  na*) 

Taap.  incraaaa  in  active  layer  for  tino  conatant  Tl,  dag.E  *) 

T«i^.  incraaaa  in  activa  layer  for  tiaa  conatant  T2,  dag.  E  *) 

TO  la  the  scale  factor  (in  dag.  E)  for  tba  shift  in  *) 
currant  threshold  for  change  in  taMparatura  *) 

aat  jen  to  1  to  include  ahift  in  currant  thraahold  due  to  taap.  ,*) 
aat  to  0  to  ignore  *) 

uaTEAJ  •) 

which  way  to  dafino  bandwidth,  1  •  loraotsian,  2  •  Modified  Loranrtlan 

gain  apactral  half  width  for  X«oraatiian  in  nicrona  *) 

gaio  apactral  halfvidth  for  aodifiad  Lorantaian,  in  aierona*) 

gain  apactral  halfvidth  for  modified  X<orantiian,  aaoond  parameter  *) 

wavelength  at  peak  of  distribution  *) 

gain  apactral  peak  offset,  in  mierona*) 

to  make  aymmatric,  make  v«0,*) 

v^altaZ«e/2  aata  oaatar  of  gain  curve  bat  ween  2  modes*) 
nuMbar  of  modes  -  2Mtl,  marim»»  is  11*) 


*) 


*  aat  u-0  to  ignora  apactral  ahift  for  carrier  daOijity  *) 


(*  n  -  electron  density  ,  t  *  time,  na  *) 

(*  CO,  Cl,  C2.  etc.  ara  the  gain  factors  for  tba  individual  nodes  *} 

(*  aO,  al,  s2  etc.  era  tba  photon  daaaitiaa  of  the  individual  modes  *) 

Snerx*  *  plaoek  e/lnmbdap  (*  anargy  par  photon,  J/photon  *) 

-19 

2.44803  10 
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ttoafp— aottftdl  (tta—nd"tijMf tart)/. 351  <*  approx,  ao.  of  xtopo  uaod  by  •laulotloa*) 

57 

doltaLe«laabdxp*2/ (2  ng  L)  (*mod«  ooporotioa,  la  microax*) 

0.000312189 

DT  ■  Toa*  (K1  (l-lxp(-t/Tl] }  ^K2  (l-Cxp(-t/T2) )  )  cbAago  ••  •  fea  of  timo  duxiag  puloo*) 

0 

DT2  -  jon*(ia(l-txpt-t/Tl])+lt2(l-lxp[-t/T21)) 

(*tKx««hoId  ohl£t‘*loi^«roturo  cbxaQo  ax  a  fuactioa  of  tiao  durlag  tha  pulaa*) 

gaaBa«a[(K  ata*Xaabdap*4)/ (4  Pl*2  aq  ar^2  dthickaaaa  width  L  BO)}  (*  apoataaaoua  amiaaioa  factor 

0.000114056 

atraaawtraaa  Kzp(DT2/TO] 

29 

2O0OOO0QOOOOOOOQQ 
b«ta2^«ta2p  Cxp{-DT2/?01 

13 

4.4  10 

Aw((bata2  ••  dthiekaaaa) /tap)  (c/ng)  (*atl«ulotad  aaiaaioa  factor, aicroaa*3/aa*) 

0.0084 

r«(r«fl  raf2)*.S  (•raflactivity  of  tha  oirrora  *} 

0.533854 

alpha (alphao4 (l/L)  Log[l/r))l  (*  loaa  coafflclant  *) 

0.00499211 


tp*ag/(e  alpha)  («  photoe  llfatlAo,  a«  •) 

0.00235038 

tp2««m/ (o  ((l/L)  ZiOgCl/r]))  (*awara9o  tlao  photoaa  apaad  la  cawity*) 

0,00560837 

jfwcorraat/ (t  width)  (*  currant  daaaity,C/aa/ialeroB*2  *) 


-14 

1.11111  10 

j*(j2/ (dthiekaaaa  aa) ) *A*tap*tp  (*ooEmaXiaad  aXactroa  danaity, aXactroaa/ (aieron^B  nx)  *) 


9‘».  9323 


ao«B( (ntraaa/aa^tap/dthickaaaa)  A  tp) 

(*traaaparancy  currant  danaity  for  deminant  aoda.  with  tanparatura  affacta  *) 
MCX-XO'^S  width/.  97  (L  traoa) 

0.0269072 

IO(2«X0*9  width/,97  (X.  aXphao/(ata  bata2)) 

0.027179 

IOO«10^9  width/. 97  (Log [ 1 /r ] / <ata  bata2)) 

0.0196074 

Ith^OQ  <IO<2  iiOO  (*cxXcuXatad  thraahoXd  currant,  aapa  *) 

0.0736936 

UP*  u<X'a/Do)  (*  ahlft  in  gain  paak  dua  to  carxiar  danxity  changa  *) 

0 

Lap[x^)  :•  UU-'h  OT  ♦x  daltxLc  (*  diffaranoa  batwaao  aoda  waaaXangth  and  gain  paak 

*(■_)  :-l£lBanatyp*<1.5, 1/ (1 ♦ (Lap [x} /BO) “2) , X/ (1* (Ltp(*) /Bl) -2* (Lwptxl /t2) ‘4) ] 

(*  danoaiinator  of  Xina  ahxpa  factor  *) 


DO  •44(0] 
1  . 


•) 


GOaDO  (n-no) 

1.  {-12.7802  •  n) 
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0.999844 

Cl-If  [ll>0,01  (D-Ba).Ol 
0.999844  (-12.7802  •  -) 
D2-If  [lOO.WI-l]  ,0] 
0.999844 

62-Zf[)O0,02  (B-BO),OJ 

0.999844  (-12.7802  •  r;] 

D3«zf  [iol.w(2),a] 

0.999377 

e3-If[iol,D3  (n-ao),0| 
0.999377  (-12.7802  •  n) 
D4-lf[K>l,*[-21,0| 
0.999377 

e4alf  (lOl ,  B4  (n-Bo)  .  0 1 

0.999377  (-12.7802  .  n) 

DS-Zf [M>2.M(3I,0) 
0.998599 

eS«ZZ(K>2,DS  (B-BO),0| 
0.998599  (-12.7802  ♦  r) 
OC-Zf  [I02,>(-31 ,0j 
0.998599 

C<aZf[l(>2,0«  (B-Be),a) 

0.998599  (-12.7802  •  nl 

B7«Zf  [103.  *(41,0] 

0.997511 

S7-Zf  [103 ,  B7  (B-BO)  .  0 1 

0.997511  (-12.7802  •  ',1 

Bt-Zf  [103,«[-4|.01 

0.997511 

Of-Zf  [103 .  Bi  (B-Bo)  ,  0 1 

0.99751  1  (-12.7802  •  r.) 

B»-iZf[lO4,N(3),0| 

0.996117 

6»-Z{[104.B9  (B-Ba),0| 
0.996117  (-12.7802  •  -I 
B10-Zf[104,N(-S|,0) 

0.996117 

el0aZZ[lO4,B10  (n-no|.a| 
0.99611 7  (-12.7802  ■  - 
Bll-Zf[10S,W(4]  .0) 
0.994418 

eil>Z2[ioS.Bll  (a-Bo).OI 

0.994418  (-12.7802  • 

B12>Zf  [10S,M[-4]  .  0) 

0.994418 

el2>Zf[105,B12  (n-no),0| 
0.994418  (-12.7802  • 
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DL3«Zf[l&<,«(71,0] 

0.992^17 

ai3-Z/[IC»<,D13  (B-DO).O] 

0.992417  (-12,7002  ♦  n) 

014«Xf  [IO<,1I[-7],0] 

C. 992417 

Gl4«Zf[l&<,D14  (n-Bo).OJ 
0.992417  (-12.7802  *  n) 

D15«If  [l&7,Wtl]  ,0] 

0.990110 

C15-If [107,0X5  {a‘Bo),01 

0.990110  (-12.7002  ♦  n) 

DlC-Zf [107, «[-•],  0] 

0.990118 

61C«Zf [107,0l<  (a-ao),0] 

0.990118  (-12.7802  •  n] 

D17«Zf  [IOt,W[9]  .0] 

0.907526 

Gl7«Zf[l&i,DX7  (B-no),0J 

0.907526  (-12.7802  *  r.) 

DXt-Zf  [IC>«,V(-9],0) 

0,987526 

Gl9«Zf[lO8,01S  (o-Bo),01 
0.907526  (-12,7802  ♦  n) 

01»-Zf(l&9,M[10].0] 

0 

Gl9«iZf[109.  019  (0-00),  0] 

0 

D20«Zf  [109.V[-10I  ,0) 

0 

G20«Zf(IO9.  020  (Q-oo),OI 

0 

D21«Zf  [10X0,11(1X1 ,01 

0 

G2X-ZC[IOlO,D2X  (n-ool.OI 

0 

D22-Zf(IOlO,ir[-lX1.01 

0 

622-Zf  [10X0,022  (n-Do).O] 

0 

(aO^»l>a2«a3^«44a5^a6-fB7««8t«9^alO^«ll«rl«r2^r34r4^r5+r64r7tr8-fr9«rl0^rll)/(t«p  A) 

at«A  tp  {1/(A  tp) ♦bo/ (A  tp)) 

13.7002 

BBOaOO 

30 
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((j*B-(60  aO^Cl  9I*  62  rl463  m2^  r2^  »2* 

0<  r3467  •446«  r446»  aS+ClO  rS+611  a44Gl2  r«4-613  *74614  r7461S  *14614  rS4 
617  *»4614  rft461»  *104620  rl04621  *114622  rll) /ff*!**) /t«p) , 

(((•t*  DO  gmmmM  a)4(*ta  GO  *0/6*10*) -aO) /tp) , 

(((•t*  01  gmamm  n)4(*ta  61  al/galaa} -«1) /tp) , < ( (at*  02  gmmmm  a)4(«ta  62  rl/6ain*)-rl)/tp) , 

(((ata  D3  6BH**  n)4(*ta  63  *2/gaiaa) -*2) /tp)  ,  <  <  (at*  D4  n)4(ata  64  r2/gain*)  ~r2) /tp) , 

<((ata  D5  gaa*M  n) 4 (ata  6S  *3/galna) •*3) /tp) , ( ( (ata  D4  gasMa  a) 4 (ata  64  r3/gaia*) -r3) /tp) , 

(((ata  D7  gas*  a)4(ata  67  a4/galaa) -*4) /tp) ,  ( ( (ata  04  gaaMa  a)4(*ta  61  r4/gaia*) 'r4) /tp)  , 

({{ata  04  gawa  a) 4 (at*  69  aS/gaia*} -*S> /tp) ,  { { (at*  010  gmmmm  a) 4 (ata  610  rS/galaa) -rS) /tp) , 

(((ata  Dll  gaasa  a)  4  (ata  611  aC/galaa)  ~*4) /tp) «  ( ( (ata  012  gas*  a)  4  (ata  612  r4/gaia*) -r€) /tp) , 

(((ata  013  gaaaa  a)  4  (ata  613  *7/gmia*) -*7) /tp) »  { ( (ata  014  giMi  a)  4  (ata  614  r7/gala*) -r7) /tp) , 

(((ata  015  gavaa  a)4(*ta  615  *1/6*1**) -*4) /tp) ,  ( ( (ata  014  gaM*  B)4(*ta  614  rl/gaia*) -rS) /tp) , 

(((at*  017  gaaau  a)4(ata  617  aO/gala*}  ~*9) /tp) ,  ( ( (ata  Dll  gaM*  a)4(ata  614  rO/gain*) -r9) /tp) , 

(((ata  019  gaaaa  a)  4  (ata  619  *10/gaia*) -*10) /tp) ,  ( ( (ata  020  gaM*  a)  4  (ata  620  rlO/gala*) -rlO) /tp) . 
(((ata  021  gasaa  a)4(*ta  621  •ll/galo*)-*ll)/tp) »  (((ata  022  ga»a  B)4(ata  622  rll/galaa) -rll) /tp)  ) , 
(t,B,a0,sl,rl.*2.r2,*3,r3,*4.r4,*5,r5,a4,r4,*7.r7,*4,r4,*9,r9,*10,rl0,*ll.rll), 

(0,a*o, 0,0«  0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 0,0),tiaaand  , 10**4, 
frograaaTraea  *>  Falsa,  MaxiauaJtapilia  *>  1., 

Zalttt*p41aa*>.01  ] ; 

(♦TlmlBg(%]  tuo  Tlaa*  3333  IM  aaooad*  •) 


fa-Ll*tFlot[tabla(  (f((k,l]),  f((k,2n).  (k,  1,  Laagtk(f )) ) , 

FlotJolaad  *>  trua. tlotftaaga*>( (tia**taft, 5) , (30, 40) ) ) 


•Craphlcs- 

(*Uisa£C(Be«tap*  1,2])  *)  (*  tba  carrlar  daaaltx  *) 

(*  0,oao,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0  •) 

JUuBa((bata2p  aa  dtklekaa**) /tap)  (e/ag)  (**tlaal*tad  «Al**loa  factor, ade»aa**3/a*4) 

0.0084 

444-4iB[f[(3  ao*tap*,i)l,{1.3.21)|*((L  width  dthiekaaaa/at*  Eaargy)  (l/tp2  10*9)  10*3/(t*p  &u) ) 

82.647 

a*tpovar«50 
tlMvtartsO 
t  iaaaad«t  laaaad 
aot  1  ek  ••net  1  ck  * 
aotlck*2«9 

5 

ddd«Tabl*[f [[5, 11)1,(11,1,21)) 

10.147134,  38.7123,  C. 0836^46,  C.C8323U,  0.0832311,  0.0019756,  0.0019756,  0.0799328,  0.0799328, 

0.0771738,  0.0737915,  C.  073791  5,  0.0698951,  0.0698951,  0.0656036,  0.0656036,  0.0610391  ,  0.061039:.  ' ,  - f  > .  : i . 
0.0563209) 

fa-Li*tPlot(Tabl*[  (f[Ck,l]),  f  [  [k,  3)  )  •  ( (X.  width  dthlckaa**/*ta  Kaargy)  (l/tp2  10*9)  10*3/ (tap  Xu)  ))  . 
(k,l,t*ngth[£l)], 

FlotJolaad  *>  Traa. Flot4aega*>( (tlaaatart.tlaaaad) , (0, a*tpewar/2) ) , 

Tick*  *>(RaAg* (t  lB*«t*rt ,  t laa*ad,  (t lMaBd''tla**tart)  /aot Ick*] ,  B*nga(0,  a*tpowar/2,  astpowar  / 2/Dot  ick  *^ ) ) ) 


-Graphics- 
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fw»U«tVlot[TafeXtt(  f  [  [k.  3])*  ( (X*  »idtth  dthickiM««/«tA  iMrgy)  (l/tp2  10^9)  10*3/ (t«p  Au) )  h 

(k,  1,  XADgth[fn  1 « klet  Joined  ->  Tru*,  tlotBAa4**><  (tlMatart,  tiM«ad/Doticks) ,  (0,  ••tpo««r/2) )  ] 


-Graphlcs- 


£n*I*l«t9lQt  [Tmbl«(  (f[[k,l]],  f[(k,4]]*((L  width  dthickn««a/*ta  Xn*rgy)  (l/tp2  10*9)  10*3/ (twp  Ajcx))), 
(k,l,l«iigth[f])], 

tlotJeiawd  ->  Tru*.  f lotRAag«~>(  (tiM«t*rt, tLM«nd) .  (0,  ••tpowwr/2) ) , 

Yieka  •>{RAag«  [timaatart.  tiaaaftd.  (tlMaad^tlMatart)  /aoticka) ,  hanga[0,  •atpo««r/2,aatpowar/2/aoticka21 ) ) 


2St 


-Graphics-* 

fn2«tlatPiotCTabl«[  (fC(k,l]l.  f  [  {k,  51  ]•(  (L  width  dthlckA*aa/ata  taargy)  (l/tp2  10*9)  10*3/ (tap  Au) )) , 

(k,l,L«ngth(fni. 

tlotJoiaad  ->  Trua.  9lotkaog«->(  ( tiMatart .  tiMwad) ,  (0.  •atpow«r/2 ) ) . 

Ticka  «>( fLanga(tiJMatart, t iaaaad,  (tiAaaad-ti*aatart) /aoticka) , kaa9a (0 , aatpowar/2. aatpowax/2/&ot icka2 ) ) ] 


-Graphics- 

fa<-ZdatriQt(Tabla[  (f[[k,l]],  f(  [k,  C)  )•(  (L  width  dthicfcaaaa/ata  Cnarvy)  (l/tp2  10*9)  10*3/ (tap  Aut  M  . 
(k,l.taBgth(fl )1. 

VlotJoiaad  ->  Trua. Flottaaga->( (tljaaatart. tiaaand) . (0. aatpowar/2 ) ) , 

Ticka  ->{kaaga  (t  l»aatart .  t  isaaad,  (tiaaaod-tiaaatart)  /aot  icka ] ,  Banga[0,  aatpowaa/2.  aatpowar/2/Bai  i  ck  a2  }  )  ) 


-Graph! cs- 

fo-Uat9lot[TablaC  (f[[k,l]1,  f  [  [k,  7]  ] «  { (X.  width  dthicknaaa/ata  laaryy)  (l/tp2  10*9)  10*3/ (tap  Au)  M  . 
lk.l,Laagth(£] II, 

PlotJoiaad  ->  Trua, PletAaBga-> ( (tiaaatart. tiaaaad) .( 0. aatpowar/2 U . 

Ticka  ->( Aanga (t i«aatart . t iaaaod.  (tiaaaad-tiaaatart ) /aot icka] , kaagafO, aatpowar/2. aatpowar/ 2/ oot i ck *2 ] ) ] 


2&t 


-Graphics- 
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Ifltk.m.  *(tk,ill*((l  En*rw)  (1/tpJ  10*»)  10*3/ (t«p  A«») ) ) , 

(k,l,t<n«th(fni, 

tlotJaiaad  ->  Tru*.  »latEuiga->(  (tiaaatart.tiMwul) .  (0.  ••tpa«*r/2)  | , 

Tlelia  t0i***t*rt ,  tiji*«ad,  (tiA*«nd-tiMatart)  /aotlcka] ,  EAng*[0.  •■tpow*r/2,  •■tpoii«r/2/notlek»21  J 1 


-Graphics- 

fwUattlscCtkblat  (ftU.lll.  f  ( (k.  t|  |  •  ( (t  •t4th  dthieko.»»/«t«  ti>«rw)  (l/tp2  10*»)  10*3/ (t.p  k»») )  I , 

(k,l,Laa«tk(fni. 

rloCJelsad  ->  Tr««,  Hot««ng»->(  .  10,  ••tpo»»»/2) ) , 

tlok»  ->(ll»ng«[tl»«»t«rt.tl»««iid.  /ootlck*) ,  »*ng»t0,  ••tp«>»«r/2,  ••tpo<»»r/2/notlck«21 1 ) 
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